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ABSTRACT 

Context. Processes driving mass assembly are expected to evolve on different timescales along cosmic time. A transition 
might happen around z ~ 1 as the cosmic star formation rate starts its decrease. 

Aims. We aim to identify the dynamical nature of galaxies in a representative sample to be able to infer and compare 
the mass assembly mechanisms across cosmic time. 

Methods. We present an analysis of the kinematics properties of 50 galaxies with redshifts 0.9 < z < 1.6 from the 
MASSIV sample observed with SINFONI/VLT with a mass range from 4.5 x 10 9 M to 1.7 x 10 11 M and a star 
formation rate from 6 Mq yr _1 to 300 M© yr _1 . This is the largest sample with 2D kinematics in this redshift range. 
We provide a classification based on kinematics as well as on close galaxy environment. 

Results. We find that a significant fraction of galaxies in our sample (29%) experience merging or have close companions 
that may be gravitationally linked. This places a lower limit on the fraction of interacting galaxies because ongoing 
mergers are probably also present but harder to identify. We find that at least 44% of the galaxies in our sample 
display ordered rotation, whereas at least 35% are non-rotating objects. All rotators except one are compatible with 
rotation-dominated (V m ax/cr > 1) systems. Non-rotating objects are mainly small objects (R e < 4 kpc). They show an 
anti-correlation of their velocity dispersion and their effective radius. These low-mass objects (log M sta r < 10.5) may be 
ongoing mergers in a transient state, galaxies with only one unresolved star- forming region, galaxies with an unstable 
gaseous phase or, less probably, spheroids. Combining our sample with other 3D-spectroscopy samples, we find that 
the local velocity dispersion of the ionized gas component decreases continuously from z ~ 3 to z = 0. The proportion 
of disks also seems to be increasing in star- forming galaxies when the redshift decreases. The number of interacting 
galaxies seems to be at a maximum at z ~ 1.2. 

Conclusions. These results draw a picture in which cold gas accretion may still be efficient at z ~ 1.2 but in which 
mergers may play a much more significant role at z ~ 1.2 than at higher redshift. From a dynamical point of view, 
the redshift range 1 < z < 2 therefore appears as a transition period in the galaxy mass assembly process. All the 
data published in this paper are publicly available at the time of the publication following this link: http: //cosmosdb. 
lambrate . inaf . it /VV DS -SINFONI 
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1. Introduction 

During the last decade the first observations of distant 
galaxies with Integral Field Unit spectrographs (IFU) have 
led to the construction of several galaxy samples from 
z ^ 0.5 to z ~ 3. These IFU observations have given new 
insights into the resolved physical properties of galaxy pop- 
ulations at various redshifts but a clear view of the tran- 
sition epoch between young, unstable and clumpy galax- 
ies and evolved and stable galaxies that form the Hubble 
Sequence is still missing. 

In the redshift range 0.5 < z < 3, large quantities of gas 
are converted into stars, producing a peak in the cosmic 
star formation at these epochs (see e.g. Hopkins & Beacom 



2006 



crete^ 



|Tresse et al.|2007 ). Large gas reservoirs are being ac 
d onto galaxies via various mechanisms: isolated events 



such as galaxy major and minor mergers (e. g.|de Ravel et ah] 
20091 IConselice et al.|2008| |Lin et al.||2008l |L6pez-Sanjuan| 



et al. 



2011) and cold gas accretio n along cosmic filaments, 

n Dekel 



a more continuous process (e.g 
fc Birnboim||2006| |Genel et al 



Keres et 
2008; 



al 

Dekel 



2005 
et al.l 



2009). 

Cosmological numerical simulations have helped to outline 
the importance that these processes may play in the build- 
up of the present-day Hubble sequence. Indeed, on the one 
hand, merging galaxies are a natural mass assembly mecha- 
nism expected in the ACDM framework where dark matter 
halos grow from hierarchical assembly. Some authors have 
shown that a pea k of major merger activity might occur 



around 1 < z < 2 ([Ryan et al.||20Q8 Conselice et al. 



Lopez-Sanjuan et aT|2009| ). On the other hand, co. 
accretion may play a signif icant role at z > 2 but 
be less usual at 



2008 
d 



Dekel et al. 2009 



z<l(e 

|2009[ ). This process see ms to be efficient f or 
larger than - 10 11 M ( |Bouche et aT||2QlQ ). 



Keres et 
halo 



gas 
might 



al. 



masses 



Recent IFU studies aim at understanding the role of 
these different processes in the precursors of local ellipticals 
and spirals. From the IMAGES sample at 0.4 < z < 0.75 



(Yang 
Rodru 



et al. 2008; Neichel et al. 2008; Puech et al. 2008 



odri gues et al.||2008| ), it has been shown that regular ro- 



tating disks are quite similar to local rotators and that 
merging is playing a significant role in galaxy mass as- 
sembly. At higher r edshift, the SINS sample at z ~ 2.2 
flGenzel et al. |2008| [Shapiro et al.||2008| |Forster Schreiber 
|et al.M2009[ |Snapiro et al.| |2009[ ICresci 
LSD/AMXZE sample at z 



20091 ICresci et al.||2009l), the 
3.3 ([Hneruc ci et al.TjgOlT) 



and the sample built by Law et al. 
are all containing many galaxies wit 
bulence. Based on numerical simulation results of unsta- 



(|2009|) at 2 < z 
itET 



2 < z < 3 
high gaseous tur- 



ble g as-rich disks (e.g. |Bournaud et al. 2007 Dekel et al. 
2009|), these authors claim that continuous cold gas ac- 



cretion along cosmic web filaments is playing a significant 
role in mass assembly at these redshifts because they find 
that half the disks are dispersion-dominated and that even 
rotation-dominated disks show a high gaseous local velocity 
dispersion (above 60 km s -1 ), uncommon fo r loc al galaxies 
dEpinat et al.| |2010|). iLehnert et aT] ( |2009[ ) and |Le Tiran 
|et al.| ( |201ip suggested that this high gaseous local velocity 
dispersion is not directly powered by shocks or Jeans insta- 
bilities due to cold gas accretion but might be related to a 
vigorous star formation at both high- and low -redshift. This 
was also supported by | Green et al.| ( |2010 ) who observed 
a few local analogs to these high-redshift galaxies that all 
have a star formation rate above 15 M Q /^r. [Gongarves et al. 
(2010) also observed galaxies with high gaseous local veloc- 



ity dispersion among a sample of z ~ 0.2 Lyman-break 
analogs that are also forming stars very efficiently. It has 
been suggested that this i ntense star formation could also 



Overzier et al 



2008]). 



be triggered by mer ging ( |Basu-Zych et al. 2007 2009a|b 



Among these various samples (IMAGES, SINS and 
LSD/ AMAZE), the redshift range 0.9 < z < 1.8, at the 
peak of the cosmic star formation history, is still poorly ex- 
plored. The MASSIV survey has been built to study this 
redshift range, which seems to correspond to the pe riod 
where the modern Hubble Sequen ce is being built (e.g. Bell 
u + - 1 "20041 iWilliams et al. 120091 . 



et 



alp 
Th« 



e first goal of this paper is to study the dynamical 
nature of 50 galaxies with 0.9 < z < 1.6 from MASSIV and 
the evolution with redshift of the fraction of galaxies in a 
given dynamical state. The second objective is to identify 
the main physical mechanisms responsible for the transi- 
tion between z ~ 2, where many disks show high intrinsic 
gaseous velocity dispersions, and z ~ 0.5 where the disks 
seem to be more "stable". The general presentation of the 
MASSIV sample selection and data acquisition strate gy is 
the subject of a companion paper ( jContini et al.|20l2| ). An 
analysis focused on the dynamical properties of rotators 



is presented in a second companion paper (Vergani et al. 
submitted) and the analysis of the spatially-resolved metal- 
licity of this firs t set of 50 MASSIV galaxies is discussed in 
an third paper ( Queyrel et al.||2012 ). 

The paper is organized as follows. In section [2j the 
MASSIV "first epoch" sample, the observations, and the 
data reduction are presented. In section |3j galaxy morphol- 
ogy and kinematics models are described. The galaxy clas- 
sification scheme is detailed in section [4] and the discussion 
of this classification is given in section |5j Appendix [A] con- 
tains detailed informations and comments on each galaxy. 

In this paper, we use the cosmological parameters Vt m = 
0.3, n A = 0.7 and H = 70 km s" 1 Mpc" 1 . 



2. Observations and data reduction 

2.1. Sample selection 

A complete description of the MAS SIV sample selection 
and general properties is presented in|Contini et al. (2012). 



Here we only present a brief summary. 

The MASSIV sample contains 84 star- forming galaxies 
in the redshift range 0.9 < z < 1.8. It was constructed from 
the VVDS (VIMOS VLT Deep Survey) spectroscopic sur- 
vey ( |Le Fevre et al.|2005[ ). This survey is I-band magnitude 
limited and compl ete up to magnitude Ia b < 24 for VVDS- 
Deep (RA= 02h, |Le Fevre et al.||2005|), I AB < 24.75 for 



VVDS-Ultra-Deep (RA= 02h, |Cassata et al.|20llj Le Fevre 
et al. in prep) a nd I AB < 22.5 for V VDS- Wide (RA= 14h 



and RA= 22h, |Garilli et al.| |2008). VVDS provides low- 
resolution (R rJ 230) spectra that lead to accurate spec- 
troscopic redshifts for 4446 galaxies in the redshift range 
0.9 < z < 2. 

The selection of the MASSIV sample was based on 
star formation activity traced by [O n] A3727 emission line 
equivalent width for galaxies with z < 1.46 and from their 
observed photometric UBVRIK spectral energy distribu- 
tion for galaxies with z > 1.46. These criteria ensure that 
the selected targets are star-forming galaxies for which 
strong emission lines can be studied using SINFONI to trace 
the kinematics. The [O n] A3727 selection criteria has been 
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tested on a pilot sample (Epi nat et al.|2QQ9 ) and has proven 
to be very efficient. Resulting from this selection function, 
the MASSIV sample provides a good representation of star- 
forming galaxies with SFR > 5 M yr " 1 at z ~ 1.5 in the 
stellar mass regime 10 9 - 10 11 M (see |Contini et al.||2Q12 
for a detailed study). 

The MASSIV "first epoch" sample discussed in this pa- 
per is a subsample of the MASSIV sample. It contains the 
50 galaxies with 0.9 < z < 1.6 (the median redshift is 1.24) 
that ware observed before January 2010. Except for two 
galaxies observed in the VVDS-Ultra-Deep, the galaxies 
studied in this paper are from the VVDS-Deep and VVDS- 
Wide samples. 



2.2. Observations 



SINFONI ( |Eisenhauer et al. (20031 |Bonnet et aL|2004| > was 
used to obtain 2D spatially resolved spectroscopy in the 
NIR of the MASSIV galaxies around the Ha line or around 
the [Om] AA4959,5007 line for four galaxies. Out of an 
initial sample of 50 galaxies no line was detected for only 
four galaxies. The success rate is much better when the 
Ra line is targeted (44/46). VVDS220148046 was observed 
at z = 2.244, whereas it was expected at z ~ 1.371 from 
the VIMOS spectrum. The redshift determination for this 
galaxy is based on the observation of [O ill] AA4959,5007 
and H/3 lines in our SINFONI data. The observations were 
obtained in service mode from period P79 to period P82 as 
part of the Large Programme 179.A-08 23 (P.I. T. Contini) 
and during a pilot program presented in |Epinat et al. ( 2009 ) 
(ESO runs 75.A-0318 and 78.A-0177). Galaxies werFob- 
served in two or three observing blocks. 

For the galaxies with z < 1.1, Ha was observed with the 
J grism over the spectral range 1.08-1.41/xm, whereas for 
galaxies with z > 1.2, Ha was observed with the H grism 
over the spectral range 1.43-1. 86/im. [Om] AA4959,5007 
was observed in the J-band except for VVDS220148046 for 
which it was observed in the H-band. SINFONI spectral 
resolution in J- and H-bands reaches ~ 2000 and ~ 2500, 
respectively. 

For seeing-limited observations, we offset the target al- 
ternatingly from one corner of the field-of-view of the in- 
strument to the opposite one (object nodding) to maximize 
the observing efficiency during the nights. This observing 
strategy allowed us to avoid sky frame acquisition. For AO 
observations, we used one sky frame for four objec t frames 
to maximiz e the time spent on sources (details in Contini 
et al. 2012[ ). In addition, we also applied a sub-dithering 
to avoid the return of the target to the same position on 
the chip. To allow for an accurate on-source pointing of 
our galaxies, we acquired them through a blind offset from 
a bright nearby star (PSF star used to measure the spa- 
tial PSF) to our target. We also observed standard (STD) 
stars for flux calibration during the same night. Individual 
exposures were 300s, 600s or 900s with a total on-source 
integration time that ranges between lh and 2h. 

Most of the data were observed in seeing-limited mode 
using the 0.125" x 0.25" pixel scale leading to a 8"x8" field 
of view with a mean seeing of 0.68 ± 0.12", considering 
only detected galaxies. However, a subset of seven galax- 
ies were observed with the laser guide star (LGS) adap- 
tive optics system using the 0.05" x 0.10" pixel scale with 
a 3.2" x 3.2" field of view. Two of them were not detected 
(VVDS020126402 and VVDS220071601) and the PSF star 



was missed for VVDS220386469. The mean spatial reso- 
lution for the four other galaxies is 0.23". The observing 
details along with the periods at which the galaxies were 
observed are listed in Table [TJ 



2.3. Data reduction 

The data reduction was p erformed using the ESO- SINFONI 
pipeline (version 2.0.0, |Modigliani et ah] |2007[ ) comple- 
mented with additional IDL and PYTHON routines to per- 
form the data processing homogeneously among the reduc- 
ers and to improve some reduction steps. 

First, the PSF stars and STD stars observations were 
reduced using standard data reduction. 

Then, for each object science frame, the following steps 
where applied: 

1. bad line removal in the raw data; 

2. dark-current, sky-background and night-sky line sub- 
traction from the raw data using the contiguous frame 
(with the target in the opposite corner due to the ob- 
serving strategy or with only sky for AO data); 

3. flat-field correction using an internal lamp; 

4. wavelength calibration using arc-lamps; 

5. flux calibration using the STD telluric star spectrum 
with standard techniques to convert counts into flux 
units. This enables one to correct for the atmospheric 
transmission and instrumental response; 

6. cube reconstruction (with and without sky lines) in 
counts; 

7. adjustment of the astrometry to match our I-band refer- 
ence images (CFHT12k survey, [McCracken et aL|2003 ) 
using the acquisition PSF star and the offsets of ob- 
ject observations with respect to this star. This method 
was not possible on AO data due to missing meta- 
infor mat ions. For some galaxies of the pilot program, 
no PSF star was observed and in other cases, PSF stars 
were saturated in the CFHT images, which provided a 
deteriorated astrometry. 

In the pipeline, the spaxels are resampled so that they are 
square (0.125" or 0.05"). 

All cubes obtained for each science frame for 
a given object were then combined using the 
SINFOUTLCUBECOMBINE recipe to obtain 
the final cubes using an additional median filtering to 
remove sky line residuals as accurately as possible from 
the final cube. 

Mono-dimensional sky spectra were extracted from the 
sky cubes. They were used to 

— determine the effective spectral resolution: using a 
Gaussian fit to approximate sky lines, the spectral res- 
olution element is found to be fairly constant over the 
whole wavelength range in the J- and H-bands and the 
dispersion of the Gaussian was estimated to <r ~ 2.8±0.2 
A; 

— quantify the noise that was considered as a Poissonian 



The data from the pilot study ( Epinat et al.|2009 ) were 
reduced again using this new procedure. 
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2.4. Map extraction 

The ionized gas kinematics of MASSIV galaxies is studied 
through the brightest emission line available in the NIR 
spectra, the Ha line or the [O ill] A5007 line in a few cases. 
IDL routines, based on the mpfit routine (|Markwardt| 



2009), were used to extract the kinematic maps from the 
SINFONI data. First, to enhance the signal-to-noise ratio 
(S/N) without degrading the spatial resolution of the data, 
a subresolution 2D spatial Gaussian smoothing (FWHM 
of two pixels) was applied on the data cubes. In addition, 
the residual cosmic rays were rejected using a 15a clipping 
on 3x3 pixel boxes. Then, for each spatial pixel, the spec- 
trum around Ha was fitted by a single Gaussian profile and 
a constant continuum (four parameters in total). To mini- 
mize the effects of noise induced by sky lines but also of sky 
line residuals on the line parameter determination, the ID 
sky spectrum was used as an estimate of the noise to weight 
the contribution of each spectral element. From these fit- 
ting techniques it was possible to recover the line flux map, 
the velocity field and the velocity dispersion map for each 
source. The instrumental spectral PSF was taken into ac- 
count to compute the velocity dispersion: <j 2 = o'lb s — o'p SF - 
During the line fitting procedure, the velocity dispersion 
was allowed to vary in the range 40 km s _1 (to avoid fit- 
ting noise) to 250 km s _1 (to avoid fitting a continuum). 
2D error maps were also derived for each quantity from the 
fitting procedure. These are statistical errors that take into 
account the error spectrum and that indicate the accuracy 
of the fit for each parameter. An S/N map was computed. 
The computation of th is map was refine d since the study 
of the pilot program ( |Epinat et al.|[2009| ). The presence of 
sky lines is now taken into account to modulate the confi- 
dence on the line detection: the inverse of the sky spectrum 
is used for weighting (w). The signal (S) is computed as 
the weighted flux of the line divided by the dispersion of 
the line (a\) multiplied by y/2ir. In the following equations, 
the weight was normalized (%2w = 1), n z is the number of 
spectral elements, AA refers to the spectral sampling and 
is expressed in the same unit as a\ : m is the modeled spec- 
trum, c is the continuum of the modeled spectrum and / is 
the observed line spectrum: 



S 



Z AA w(m — c) 



Cr A V27T 



The noise (N) is computed as the square root of the 
weighted variance of the residual spectrum: 



(1) 



This yields the S/N: 

S n z A\^w(m — c) 

N ~ ax^ir^wil-m) 2 ' 

If the weight is constant, then, according to equation [T] 
the S/N is simply the ratio of the intensity of a Gaussian 
line over the RMS of the spectrum. This criterion is ro- 
bust and an S/N threshold of three was used to clean the 
kinematic maps (see Appendix [bI). Extra- cleaning was per- 
formed manually to remove isolated pixels (less than about 
five adjacent pixels) or groups of pixels associated with data 
artifacts that could have passed the criteria. The existence 
of any I-band counterpart was checked when there were 
groups of more than four pixels. 



3. Modeling 

3.1. Morphology modeling 

The stellar continuum is barely detected in the SINFONI 
data. To study the morphology of the stellar component, 
we used the best CFHT I-ban d images available, i.e. f rom 
the CFH12K/CFHT survey flMcCracken et al.|[2003[ ) for 
the galaxies in the 14h field and from the CFHT Legacy 
Surve}Q with the best seeing fo r the galaxies in the 02h 
and 22h fields. We ran GALFI T ([Peng et al.|2QQ2| ) on those 
images using a Sersic profile ( |Sersic||1968| ) to recover the 
morphological parameters: the center, the position angle of 
the major axis (PA), the axis ratio (b/a), the effective ra- 
dius (R e ), the Sersic index (n) and the total magnitude. 
GALFIT convolves the model to the spatial PSF to con- 
verge into a set of beam-smearing corrected parameters. 

To obtain robust estimates of the parameters and of 
their associated error bars, it is mandatory to control the 
PSF as well as possible. Indeed, for objects with sizes com- 
parable to the resolution, the use of an overestimated PSF 
leads to low axis-ratios with small error bars since no strong 
elongation may be reproduced: this would lead to system- 
atic underestimating of structures in the galaxies. We are 
indeed in most cases in this situation when the size on the 
objects is of the same order as the spatial resolution for 
seeing-limited ground-based imaging surveys. Thus, in each 
CFHT field, we randomly selected stars to characterize and 
follow the variation of the PSF. The detection of stars is 
based on color and morphological criteria for the 02h and 
22h fields and on morphological criteria only for the 14h 
field. From these stars, the width, the axis ratio and the ori- 
entation of the PSF are found to vary much across the 14h 
and 22h fields. The most affected field is the 14h one where 
the seeing smoothly varies from 0.6" to 1.0". Therefore, we 
selected specific PSFs in each field, located at 45" from each 
galaxy on average. The study of the PSF distribution over 
the various fields enabled us to determine the uncertainty 
on the PSF FWHM to be - 0.1 pixel (~ 20 mas). It is 
estimated as the deviation from large-scale variations. 

GALFIT produces residual maps that where used to 
check the convergence of the fits. It was sometimes nec- 
essary to fit secondary objects to have a correct fit (cf. 
Appendi x [A| . The parameters of the fits are summarized 
in Table EP 

We tried in a second step to fit the morphology while 
fixing the position angle of the major axis to the value de- 
termined from the kinematic s modeling (cf. Sect ion 3.2) as 
was done in the pilot study ( |Epinat et al.| |2009). However, 
this leads to a bias toward high axis ratios because GALFIT 
cannot match the elongation when the position angle is 
fixed and finds that round morphologies match the data 
better. Therefore, we adopt the morphologies as derived in 
an unconstrained way from the CFHT images in the fol- 
lowing. 



3.2. Kinematics modeling 

Among the various dynamical states of galaxies, the easiest 
to probe is that of the rotating disk. We therefore tested 
the likelihood of this hypothesis for the galaxies in our sam- 
ple and recovered the fundamental dynamical parameters 
within this hypothesis. The velocity field is accordingly fit- 



1 http : //www . cf ht . hawaii . edu/Science/CFHLS/ 
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ted with a model that assumes that the ionized gas is lo- 
cated in an infinitely thin rotating disk, as in Epi nat et al.| 
(2009). The rotation curve is described by a linear slope in 
the inner parts and a plateau in the outer parts. The veloc- 
ity along the line of sight is computed taking into account 
geometrical position effects. The model parameters are 

— x c , y c : the center coordinates; 

— z: the redshift corresponding to the systemic velocity; 

— i: the inclination of the gaseous disk; 

— PAk'. the position angle of the major axis; 

— Vt'. the plateau rotation velocity; 

— r t : the turnover radius at which the plateau is reached. 
The method used to adjust the models is described in 



detail in Epinat et al. (2010). It is based on a \ mm 



imization and takes the velocity error map into account 
to minimize the contribution of the regions with low S/N. 
The spatial PSF is taken into account in these models and 
is described with a 2D Gaussian. Its FWHM is computed 
on the PSF stars associated to each observation. To com- 
pute the model velocity field, a higher resolution velocity 
field is constructed (with at least eight pixels in the PSF 
FWHM) from the analytical model and a high-resolution 
line flux map has also to be built. Indeed, in a final low 
spatial resolution element, the contribution of the line of 
sight velocity at higher resolution is weighted by the true 
line flux distribution. We used a linear interpolation of the 
observed flux map to avoid making any assumption on the 
real flux distribution. This is one of the major uncertainties 
of our models because the real line flux distribution could 
be more clumpy than observed. 

Our model allows us to compute the circular velocity 
within the disk hypothesis. This assumption is not realis- 
tic for all objects in the MASSIV sample in which we also 
expect mergers (ongoing or late stage), spheroids, or struc- 
tures with chaotic motions. However, this hypothesis allows 
us to compute a map that contains only the beam-smearing 
effect on the velocity dispersion, due to the blurring of 
large-scale motions. Thus, by subtracting quadratically the 
map deduced from the model to the observed velocity dis- 
persion map, a velocity dispersion map corrected for the 
beam-smearing effect is obtained (noted "a residuals" in 
Appendix [B]) . This correction is also valid at first order for 
non-rotating objects, even if it can overestimate the correc- 
tion in the inn er parts. Details of t he method are given in 
Appendix A of Epinat et al.|(|2010 ). The maps of the mod- 
els are shown in Appendix |B| and the resulting parameters 
are given in Table p] For non-isolated galaxies (see Section 
[4|, we fitted the various components separately when they 
were sufficiently extended (VVDS220397579). In the table, 
companions have the suffix "s". 

The model rotation curve only reproduces the velocity 
field and is not based on a gravitational potential model (as 
for instance in |Forster Schreiber et al.||2006| |Cresci et al.| 



2009 Gnerucci et al.||20TT) since the shape of the gravita- 



tional potenHaTisiinknown and can probably not be de- 
scribed by a stellar component only. Note that a gaseous 
thin disk in rotation is not incompatible with a spheroidal 
stell ar distribution su ch as was observed in local ellipticals 
(e.g. Sarzi et aL|2006 ) or as suggested from numeri cal sim- 
ulations (e.g. |Bournaud et al.||2007| |Xu et al.||2010 ) . 

As described in the following sections, some parameters 
are difficult to constrain from the kinematics. To reduce 
the number of free parameters, we constrained the center 



and the inclination from the morphology to model the kine- 
matics, assuming that the stars and the ionized gas follow 
a common distribution. Thus only four parameters remain 
free and can be reasonably constraine d from our observed 
velocity fields. Using these constraints, Epinat et al. (2010) 
have shown from 137 galaxies of the GHASP sample (the 
largest 2D k inematics sample of nearby late- type galaxies 
described in |Epinat et al.| |2008b a) projected at z ~ 1.7, 
that this method enables one to recover the other model 
parameters statistically. 

3.2.1. Center 

Owing to the low spatial resolution of our data, the kine- 
matic center of a given galaxy is b arely constrained from 



the kinematics ( |Epinat et al. 20101). Therefore the centers 
are taken from trie I- band morphology since we expect that 
in the inner parts of a galaxy at these redshifts, the stars 
dominate the gravitational potential. 

Thanks to the method we used, the astrometry in 
the SINFONI data cubes matches the I-band astrometry. 
However, a post correction was applied to match the outer 
isophotes of the galaxy in the CFHT images and in the Ha 
maps. The median offset of this post correction is 0.18", 
which agrees with the accuracy of the SINFONI pointing 
system, which is estimated to be 0.1-0.2" (see SINFONI 
manual, using a guiding star instead of offsetting from a 
bright star). This is also the final accuracy of our astrom- 
etry. This offset was computed for 38 galaxies. The other 
galaxies were excluded (i) when no PSF star was observed 
or (ii) when they were observed with AO because the ob- 
serving sequence did not allow us to compute an astrometry 
correction. We find a good agreement between the centers 
derived from the I-band images and the peak in the Ha 
maps. 



3.2.2. Inclination 

Disk inclination is a critical parameter to estimate because 
it is directly linked to the rotational velocity (V#) of the 
disks. Indeed, we measured the velocities projected along 
the line of sight: V\ os — Vq x sini Consequently, the ob- 
served velocity has to be corrected for the inclination and 
this correction is larger for galaxies with low inclination. 
Owing to this degeneracy between inclination and rotation 
velocity in rotating disk models, the inclination is not well 
constrained from the kinematics alone. This degeneracy can 
theoretically be solved for high-resolution observati ons but 
not for ob servations with strong beam-smearing (Epinat 
et al.pOlO ): this would lead to an erroneous rotational ve- 
locity. Therefore the inclination was constrained from the 
morphology axis ratio. In addition, since the uncertainty on 
the inclination (its sine) directly impacts the uncertainty on 
the rotation velocity, we paid special attention in deriving 
realistic error bars on the inclination. The thickness of the 
disk, which is considered as null in our models, could im- 
pact the determination of the inclination, mainly for edge- 
on galaxies. However, for these objects, the deprojection 
has a negligible impact on the determination of V max (less 
than 3% considering a thickness leading to b/a ~ 0.2). 

To take into account both the uncertainty on the model 
and the uncertainty on the PSF used to recover galaxy pa- 
rameters in GALFIT, we used a Monte Carlo method. For 
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Fig. 1. Example of the method used to estimate the uncer- 
tainty on the inclination for galaxy VVDS020164388. The 
x-axis corresponds to the input inclination of the modeled 
galaxies, the y-axis corresponds to the inclination recovered 
by GALFIT for each modeled galaxy and the red dashed 
line corresponds to the inclination of the real galaxy. The 
uncertainty is measured along this line as the range that 
contains 68% of the points (lcr). 



each galaxy we simulated 2000 synthetic sources with the 
same structural parameters (magnitude, Sersic index, effec- 
tive radius, position angle), but a random value of b/a. The 
PSF to simulate these sources was also randomly chosen 
among different stars in the fields. Poisson noise was added 
and the simulated sources were placed in a nearby piece of 
sky from the original image (three pieces for each galaxy). 
The simulated sources were then modeled with GALFIT us- 
ing the PSF used for the real galaxy. The final uncertainty 
was estimated from the range of input inclinations that led 
to the inclination that was measured with GALFIT on the 
real galaxy (see Figure [T]). For the smallest galaxies, the in- 
clination was not constrained enough. We therefore decided 
to use an inclination of 60° , the median value for randomly 
distributed disks, and an uncertainty of 24°, which yields a 
probability of 0.68 (la). 

We also checked that there was no bias by study- 
ing the inclination distribution for the MASSIV sample. 
Theoretically, for uniformly randomly oriented thin disks 
(i.e. with a null thickness) the probability to observe a 
disk with an inclination between 0i and 62 is equal to 
I cos#i — cos O2 |, leading to distribution as displayed in 
black in Figure [2] with a median value of 60°. In Figure [2j we 
observe that the distribution of MASSIV galaxies with con- 
strained inclinations (filled blue histogram) misses face-on 
objects (i = 0°) and edge-on objects (i = 90°) and has an 
excess of objects with intermediate inclinations compared 
to the theoretical distribution. The most plausible expla- 
nation for the lack of edge-on morphologies is that disks 
are thick. It could also be that extinction is higher in these 
galaxies, inducing a non-detection of the [Oil] A3727 line 
in the VVDS spectra. Figure [l] clearly shows that GALFIT 
cannot recover extreme inclinations. However, the lack of 
face-on objects could also be attributed to the small num- 
bers expected or to a computing bias. It could also be 
that galaxies do not have regular morphologies and that, 
due to surface brightness dimming, one can only observe 
the clumpy irregular emission. Indeed, I-band morphology 
could be in some cases contaminated by gaseous emission 




20 40 60 

Inclination (degree) 

Fig. 2. Distribution of the inclinations for the 40 galaxies 
of the MASSIV sample (with a constrained inclination) us- 
ing a Sersic profile with GALFIT and assuming they are 
thin disks (filled blue histogram) compared with the theo- 
retical distribution for a randomly oriented thin disk (black 
histogram) . 



lines like [Oil] A3727 because the MegaCam/MegaPrime 
i' band filter covers the 0.70 to 0.84 jim spectral range, 
which corresponds to [0 11] A3727 redshifted between 0.87 
and 1.26. Unfortunately, the CFHT I-band imaging res- 
olution does not allow us to conclude about the clumpy 
emission. The peak around 60° could also be attributed to 
(i) the thickness of disks and (ii) the observation of galax- 
ies that are not disks (mergers for example). The median 
of the distribution for our sample (59°) is compatible with 
the median expected for a distribution of disks with random 
inclination. 



3.2.3. Constraints on the other model parameters 

Since inner velocity gradients reaching a plateau (V t ) within 
less than one pixel could not be resolved with our spatial 
resolution which is higher than four pixels, the turnover ra- 
dius r t was constrained to be at least one pixel during the 
fitting process. However, when the model converged toward 
this value, the statistical error becomes null. On the other 
hand, when the plateau is not reached according to the % 2 
minimization, the error on r t becomes large and so does the 
error on Vt because the models mainly constrain the slope. 
Therefore the errors on these parameters are difficult to in- 
terpret and to propagate to compute the error on the inner 
slope of the velocity gradient. A solution would consist in 
using a model described by the slope a and the turnover 
radius r t to estimate the error on the slope. This test has 
been performed but did not give more realistic errors in 
many cases. Because this solution was not convincing, we 
used the original model. For these galaxies, no error is in- 
dicated for r t and V t in Tabl e [3} Although r t is not well 
constrained for small galaxies, Epinat et al. (2010) showed 
with 137 galaxies projected at z ~ 1.7 that leaving r t as 
a free parameter statistically gives a good estimate of the 
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shape of the rotation curve. It is therefore necessary to let 
r t be a free parameter, in particular for the largest galaxies. 

3.2.4. Parameters deduced from the kinematic maps and 
models 

Radius of the Ha extent Ri as t 

This radius is computed from the cleaned maps. The center 
and position angle from the kinematic best-fit models were 
used to derive the radius of each pixel. Ri as t is the radius 
that both sides of the galaxy reach. 

Velocity shear V shear 

V s hear was computed as the total shear observed in the 
modeled velocity field shown in Appendix [B] (i.e. with 
an S/N threshold of 3) which takes into account the 
uncertainty on the observed velocity field. It is a projected 
velocity along the line of sight and is not inclination- 
corrected and consequently makes no assumption on the 
geometry of the galaxy. The use of the model enables us to 
smooth the velocity field at the outskirts. 

Maximum rotation velocity V max 

This was computed according to the model at the Ri as t 
radius. Two sources of uncertainties were added in 
quadrature to compute the final uncertainty on V max . 

The first one is the uncertainty on the inclination (cf. 
Section 3.2.2). To propagate uncertainties from inclination 



to rotation velocity, we used a Monte Carlo method assum- 
ing a Gaussian distribution for the inclination. The uncer- 
tainty on Vmax was the standard deviation of the resulting 
distribution on this parameter. 

The second source of uncertainty is related to the mod- 
eling. Because the uncertainty on Vt is a statistical one, it 
can be fairly small. Instead of usin g an approach base d on 
these errors, the GHASP sample (Epinat et al. 2QQ8a|b ) was 
exploited to compute model uncertainties related to the size 
of the gala xies with respect to the seeing. Using the same 
method as |Epinat et al.| ( |2Q1Q| ), 136 GHASP galaxies were 
projected at z = 1.33 under typical seeing and sampling of 
SINFONI observations of the MASSIV sample and the S/N 
of the simulations was adjusted to match the Ha fluxes from 
the MASSIV sample. Figure 13 of |Epinat etaL] fl2010| ) was 
reproduced with these new simulations. This figure displays 
the relative error on the maximum velocity determination 
with respect to the beam-smearing parameter defined as 
B = D25/2.S, where D25 is the optical diameter and s is the 
seeing FWHM. A linear regression was performed to esti- 
mate the evolution of the accuracy of the fit with galaxy 
size. The best fit gives 



SV m 



Vrr, 



27.5 - 5.8B 
100 



(2) 



The percentage of galaxies below this relation is about 
60%. Assuming a Gaussian distribution, 1-<j corresponds 
to a confidence level of 68%. We therefore simply ap- 
proximated that SV max is the 1-a uncertainty. For the 
MASSIV sample, B was estimated from half-light radii 
(Re ) determined by GALFIT on I-band images (see Section 
3.1). Indeed, assuming an exponential disk, the following 
relation can be written: D25/2 = 1.9R e . The model uncer- 
tainty was consequently determined following Equation |2j 
A minimum uncertainty of 10 km/s was imposed, however. 



This approach is well-adapted for rotators because it was 
built from a control sample of rotators. 

Local velocity dispersion a 

The v elocity dispersion was computed as in Epi nat et al.| 
(2009): this is the average of the velocity dispersion map 



corrected for beam-smearing effects (see Section 3.2) and 
spectral PSF. A weight proportional to the inverse of the 
uncertainty on the velocity dispersion (estimated before 
correcting for beam-smearing and spectral resolution) 
was attributed to each pixel to compute the average. 
The uncertainty on this parameter was computed as the 
weighted standard deviation. 

4. Kinematics classification 

To distinguish between galaxy formation scenarii, it is nec- 
essary to know the dynamical state of galaxies for large 
samples and at various redshifts. 

4.1. Previous classification schemes of large samples 

The first kinematics classification of distant galaxies ob- 
ser ved using inte g ral fi eld unit techniques was provided 



by |Flores et 
sample ( 0.4 



al. 



< z 



([2006} in the frame of the IMAGES 
< 0.75) w hich contains 68 classified 
galaxies ( jNeichel et al. 2008). Their classification is a 

HST optical images and 



visual one that relies on bot 



FLAMES/GIRAFFE data. The galaxies were split into 
three classes: (i) rotating disks whose morphological and 
kinematic position angles match well and with a velocity 
dispersion peak in the center, (ii) perturbed rotators whose 
morphological and kinematic position angles match well but 
that have a peak in the velocity dispersion offset from the 
center and (iii) complex kinematics systems with both a 
disagreement between morphological and kinematic posi- 
tion angles and a peak in the velocity dispersion offset from 
the center. The IMAGES survey has a majority of galax- 
ies with complex kinematics (44%), and a nearly equivalent 
fraction of rotating disks (29.5%) and of pertur bed rotators 



(26.5%) (Yang et al.||2008 Neichel et al 



2008). From these 



numbers, they concluded that at that epoch merging is still 
more active than in the local Universe. 

The classification of the SINS sample (z ~ 2.2) relies on 
the asymmetries measured on both the velocity and veloc- 
ity dispersion fields to distinguish rotator-like from merger- 
like galaxies. These a symmetries were der ived either from 
a kinemetry analysi s ([Shapiro et al . 2008) or from a qual- 
itative assessment ( jForster Schreiber et al. 2009| ). Then, 
comparing the mean local velocity dispersion and the rota- 
tional velocity, these galaxies (both mergers and rotators) 
were classified either as rotation-dominated (V max /cro > 1) 
or dispersion-dominated (V max /<Jo < 1). For galaxies with 
small size or low S/N, they instead compared the full veloc- 
ity shear v oos with the integrated line width (j int a nd used 
a threshold of v b s / (^ int) ~ 0.4. On the one hand, |F6rster 



|Schreiber et al. ( |2009 ) found that one third of the 62 galax- 
ies in the Ha SINS sample contains mergers. On the the 
other hand, they also found that one third of this sample 
contains rotation-dominated systems, another third corre- 
sponds to dispersion-dominated systems, and the last third 
are not classified. The authors interpret the large amount 
of dispersion-dominated disks as an evidence for cold gas 
accretion along cosmic web filaments. 
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For the LSD/ AMAZE sample, the classification only 
distinguishes rotator-like from perturbed galaxies because 
galaxies have a fairly low S/N. This classification relies on 
the modeling of the v elocity field as an inclined X-Y plane 



( Gnerucci et al.pOll ). The criterion is based on the x 2 of 



the fit by the plane and on a constraint on the reliability of 
the inclination of the plane. 

In the pilot run study of MASSIV flEpinat et aL|2009D , 
which was based on a visual kinematic classification, nine 
galaxies with 1.2 < z < 1.6 were classified into three groups: 
(i) three galaxies are mergers, showing disturbed kinemat- 
ics and possibly several components, (ii) two objects are 
classified as rotation-dominated disks, with clear signs of 
rotation and with a maximum rotation velocity higher than 
the mean local velocity dispersion and (iii) four galaxies are 
dispersion-dominated disks, with clear signs of rotation but 
with a maximum rotation velocity lower than the mean lo- 
cal velocity dispersion. It was concluded that about one 
third of these galaxies are observed during some merging 
event, whereas one third are rotation-dominated and the 
other third are dispersion-dominated. 

In this paper, we present a new kinematic classification 
for the MASSIV sample to put constraints on the dynamical 
state of galaxies at z ~ 1.2. The classification scheme was 
developed in several steps. First, eight people of the collab- 
oration independently defined their own criteria. These cri- 
teria were then reconciled and discussed during a common 
session and led to a unique classification based both on the 
close environment and on the velocity shear strength. The 
visual classification helped in defining measurable criteria 
to build an automatic, thus reproducible, classification (this 
led to minor changes in the final classes) . The final classifi- 
cation was also refined to describe the dynamical state and 
support of the galaxies. 

4.2. Galaxy small-scale environment 

The study of galaxy close environment is crucial for in- 
ferring a merger rate on the population probed by the 
MASSIV sample at z ~ 1.2. Using SINFONI data and 
CFHT imaging, galaxies could be classified as interacting 
or isolated. Flags have also been attributed to this classifi- 
cation to qualify its reliability. Flags "A", "B" and "C" mean 
secure classification (> 90% probability), confident (~ 75% 
probability) classification and poor reliability (~ 50% prob- 
ability), respectively. 

The SINFONI field-of-view was explored to detect emis- 
sion lines that might be attributed to companions. The 
same field was also explored systematically in the I-band 
images. Owing to the nodding strategy, the field-of-view 
is larger for seeing-limited observations. In that case, the 
shape of the field results from the superposition of two 8" 
square fields overlapped on their opposite corners (4" quad- 
rants). Thus, the field-of-view of the combined data around 
each target can be as large as 6" (corresponding to ~ 50 
kpc at z ~ 1.2). In the non-overlapping regions, the expo- 
sure time is half the total on-source exposure time but is 
sufficient to at least detect emission lines of objects. The 
explored field is not a full 12" x 12" square field but only 
78% of this area. For AO observations, the field of view is 
restricted to 3" (corresponding to ~ 25 kpc), which is insuf- 
ficient for studying the environment and therefore only the 
I-band image was explored to see if a possible companion 
was observed in the 6" around the main object. 



Galaxies were classified as interacting with a confidence 
flag "A" when the following conditions were met: (i) the av- 
erage S/N (see Table |4| has to be higher than three in an 
area larger than the seeing, (ii) the systemic velocity differ- 
ence between the two components has to be lower than 1000 
km s _1 and (iii) an optical counterpart has to be observed 
within the CFHT I-band image. If the companion was ex- 
tended but slightly smaller than the seeing, we assigned a 
flag "B". When no counterpart was observed in I-band or if 
some galaxies were observed in the neighborhood in I-band 
but not in the SINFONI datacube, the object was consid- 
ered as isolated with a flag "B". A flag "C" was assigned 
for AO observations unless absolutely no other galaxy was 
observed in I-band around the object. This is detailed in 
Appendix [Aj 

In some cases, the morphology is elongated and several 
Ha blobs can be detected along the elongation, sometimes 
also in I-band images (see comments in Appendix |A|). It 
is not straightforward to asses if these blobs are small ob- 
jects that are accreted and not star forming regions within 
one unique galaxy. A monotonic velocity field along the 
direction of the major axis is likely evidence for rotation. 
However, if the elongation is the result of the projection of 
the two components and if the two objects are small (no 
gradient detectable in each component), this could mimic 
a monotonic velocity field caused by the difference of radial 
velocity between the components. Some additional features 
can also support a merger hypothesis: (i) the angular sep- 
aration is larger than 3" (ie ~ 25 kpc), because galaxies 
this large are not expected at high redshift, nor are they 
observed in the local Universe; (ii) the radial velocity gap 
is larger than ~ 600 km s _1 , because this rotational veloc- 
ity would imply an unrealistic dynamical mass for a unique 
object. If one of these features was seen, the galaxies were 
classified as interacting but with a flag "B". If none of them 
was observed, the galaxies were classified as isolated but a 
flag "B" was assigned. For the difficult cases, we were able to 
use the morphology traced by old stars (mainly observable 
in near- infrared bands): the old star population may not 
follow the star forming regions in case of unique objects. 
The S/N of these J- or H-band continuum morphologies 
recovered from SINFONI data is fairly low, however. It is 
therefore not straightforward to arrive at definite conclu- 
sions. 

When components are very close, i.e. about to merge, 
the situation is more complex. However, perturbed line pro- 
files (e.g. VVDS020294045, cf. Appendix |A| in the region 
between the two blobs favor a strong discontinuity in the 
velocity field, which points toward a two-component hy- 
pothesis. These cases received a flag "B". 

For ongoing mergers one expects to observe peculiar 
kinematic signatures. Therefore, we classified some galaxies 
as interacting because their velocity fields, velocity disper- 
sion maps and morphologies were presenting perturbations 
in the same positions, suggesting for instance the presence 
of tidal arms (e.g. VVDS020283083, cf. Appendix 0. Since 
these signatures are more subjective, these galaxies usually 
have a flag "B" for the isolation criterion. 

Other galaxies were classified as isolated. However, 
when there were some perturbations in the velocity field 
or a peak in the velocity dispersion map, this could be in- 
terpreted as signs for merger remnants, but because it con- 
cerns the galaxy outskirts, it is more likely noise and we 
assigned a flag "B". 
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Fig. 3. Disagreement between morphological and kine- 
matic position angles (within the errors) as a function 
of the velocity field residuals normalized by the velocity 
shear. Blue dots and red squares represent rotating and 
non-rotating galaxies. Symbols with black contours mark 
interacting galaxies. The symbol size is related to the S/N 
of the SINFONI data: small, medium and large symbols 
refer to galaxies with S/N < 5, 5 < S/N < 10 and 
S/N > 10, respectively. One galaxy is not in the expected 
region (VVDS020294045, cf. Appendix [£} . 



4.3. Velocity shear strength 

For the kinematic classification, we used a first criterion 
as simple and objective a s poss ible based on the velocity 



shear V s hear (see Section 3.2.4). The sample was accord- 



ingly divided into low- velocity shear galaxies {Vshear < 100 
km s _1 ) and high- velocity shear galaxies (Vshear > 100 
km s _1 ). This gives a rough idea about the dynamical 
state of a galaxy without taking into account more com- 
plex motions. For VVDS220397579 and VVDS220544394 it 
was possible to study the velocity shear of the companion 
as well. 

This criterion has the advantage that it can be easily 
measured. However, the value of Vshear for a given galaxy 
may vary with the orientation. This is particularly true for 
galaxies in ordered rotation (V r shear is lower when observed 
face-on). However, the distribution of the inclinations, as 
seen in Figure [2j shows that only eight galaxies according 
to the theoretical distribution and only three according to 
the observed distribution may need a correction larger than 
a factor of two (corresponding to an inclination of 30°) for 
the measurement of the velocity shear. 

In galaxies with a low- velocity shear there might be sev- 
eral classes of objects: (i) very low-mass objects, (ii) face-on 
rotating galaxies, (iii) ongoing mergers in a transient state, 
(iv) spheroids, if one expects that the gas in these objects 
follows the distribution of stars, and (v) galaxies with a 
non/slowly rotating gaseous component. 

4.4. Rotating and n on- rotating galaxies 

For a specifi c analysis (e.g. Tully-Fisher relation, see 
|Vergani et al.| submitted), it is essential to be able to define 
a robust sample of rotating galaxies. We expect a rough 
agreement of the morphological and kinematic position an- 
gles for rotators. They can be slightly different due to per- 



turbations like bars or strong spiral arms. We also expect for 
these galaxies that the rotation motions dominate over per- 
turbations. To distinguish rotators from non-rotating galax- 
ies, we made a diagram in Figure [3] with two quantities that 
quantify these two arguments: the disagreement between 
morphological and kinematic position angles derived from 
the models presented in sections 3.1 and 3.2 within the 
uncertainties, \PA m - PA k \ - y/SPA^ + 5PA\, and the 
mean weighted velocity field residuals normalized by the 
velocity shear. Rotators are therefore defined as galaxies 
with a position angle difference lower than 20° and velocity 
field perturbations lower than 20%. These thresholds were 
chosen from the distribution seen in Figure [3] because they 
enable us to isolate a clear cloud of points near the origin 
(the rotators). We note that these values are also represen- 
tative of relatively small deviations from rotational motion 
and position angles. There is one exception in this diagram: 
VVDS020294045 was classified visually as non-rotating al- 
though it lies in the rotator's region of the diagram. Indeed, 
the velocity shear visible in its velocity field is probably 
caused by a very close companion that thus mimics a ro- 
tating disk velocity field. The shape of the line suggests a 
merger. If we were to exclude the companion, which is over- 
lapping with the main galaxy, the resulting velocity field 
would probably be classified as non-rotating (see Appendix 

We checked the agreement between this classification 
and the classification b ased on kinemetry used in the SINS 
sample (Shapiro et al. 2008). The agreement is reasonable 
but we find that galaxies that would be classified as merg- 
ers according to their criteria are mainly those that we 
have classified as non -rotating. Indeed, the classification of 
Shapiro et al. (2008) relies on the position in a diagram 
in which galaxies are placed according to their asymmetric 
velocity and velocity dispersion components. Galaxies show 
a strong correlation between these two quantities. This cor- 
relation arises because asymmetries are normalized by the 
global velocity shear. We are able to show that the nor- 
malized residuals of the model are well-correlated with the 
normalized asymmetries in the velocity field. This is true 
for large and small galaxies. Thus, our classification adds 
the agreement between the position angles in morphology 
and kinematics as a new criterion. One advantage of this 
criterion is that it is not restricted to galaxies observed with 
a good S/N and with many resolution elements. 



4.5. Dynamical support 

The last criterion is related to the dynamical support of 
galaxies. Indeed, it has been shown that at high redshift 
(z > 2), galaxies are on average more dispersion-dominated 
than in the intermediate (z ~ 0.5) and local Universe. One 
of our goals is to determine the evolution of the fraction of 
dispersion-dominated systems with redshift since MASSIV 
might probe a transitional redshift range. 

To quantify the dynamical support, we assumed that all 
galaxies have some rotation and that this rotation is ade- 
quately described by the rotating-disk model. Then, we ad- 
ditionally assumed that the local velocity dispersion of the 
gas is representative of the random motions in the galaxy. 
Thus, the dynamical support is estimated by the ratio of 
the maximum rotation velocity over the local velocity dis- 
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persion V m ax/&- This description is intrinsically better jus- 
tified for rotators, however. 

The various classes and related parameters are summa- 
rized in Table [U 

5. Discussion 

Using the classification scheme described above, we were 
able to classify as rotating or non-rotating the 46 galax- 
ies (including one companion) in the redshift range 0.9 < 
z < 1.6 for which some emission line was detected in the 
SINFONI datacubes. For part of the discussion below we 
limit the sample to th e 36 galaxies with an average S/N (de- 
fined in Section 2.4 and reported in Table |4| higher than 
5 for which the kinematic classification is more robust (we 
note that these low S/N galaxies are the majority of galax- 
ies with Ri as t / seeing < 1, which means that the kinematic 
model might be less robust). Accordingly among an initial 
sample of 46 secure galaxies with 0.9 < z < 1.6, 22% are 
not robustly classified or not classified from their dynam- 
ical properties. Th is proportion is on th e same order for 
the LSD/ AMAZE flGnerucci et !Ip0lT| ) . In the MASSIV 
sample, these galaxies have on average lower stellar masses 
and star formation rates (see Table pi, the latter explain- 
ing the lower S/N of these observations. We also emphasize 
that the high dispersions are caused by the two galaxies ob- 
served with AO, which have larger stellar masses and SFR 
than the rest of the unclassified galaxies. For those two 
galaxies the low S/N is due to the small pixel size rather 
than an intrinsically low SFR. Concerning the close envi- 
ronment criterion, the 41 systems classified with flags "A" 
and "B" are considered. 



5.1. Rotating disks vs. non rotating galaxies 

Rotating disks represent at least 44% (20 out of 46 detected 
galaxies) of the MASSIV "first epoch" sample and about 
56% of the high S/N sample. We find a lower percentage of 
secure non-rotating systems (35% or about 44% of the high 
S/N sample) that have no observed rotation in the gaseous 
component or have very disturbed kinematics with respect 
to their broad band morphology. The physical properties 
of these two classes of galaxies are on avera ge different. 
We have performed Kolmogorov-Smirnov tests ( |Press et al.| 
1992 ) on stellar mass, SFR, half-light radius and velocity 
dispersion distributions. The two populations look different 
in terms of stellar mass and SFR at a level slightly higher 
than la (the probability to follow the same distribution is 
lower than 0.2). Distributions for size and velocity disper- 
sion are fully compatible. Rotating objects are on average 
more massive (M star = 4.0 x 1O 1O M ), more star- forming 
(SFR = 6OM yr- 1 ) and have larger radii (R e = 3.8 
kpc) than non-rotating ones (M star = 1.6 x 1O 1O M , 
SFR = 39M yr" 1 and R e = 2.7 kpc). On the other hand, 
these two types of galaxies have very similar velocity disper- 



sions around 60 km s 1 (see Table [5|. At z ~ 2.2, Forster 



Schreiber et al. ( 2009 ) already noticed a similar trend from 
the SINS sample: dispersion-dominated systems are on av- 
erage smaller than rotation-dominated ones. 

The percentage of rotating systems is higher in MASSIV 
1.2 than in the LSD/A MAZE sample at z ~ 3.3. 



at z 



Indeed, |Gnerucci et al. ( 2011[ ) found a lower limit of 34% 
of rotators in their sample (11 out of 32 detected galax- 
ies that they were able to classify). Their classification is 



closest to ours since it relies both on velocity field mod- 
eling and on the agreement between morphological and 
kinematic position angle of the major axis. However, owing 
to their small statistics and their different selection func- 
tion it cannot be excluded that these proportions are com- 
patible. At z ~ 2.2, the comparison with the SINS sam- 
ple is more difficult because of the selection function, as 
for LSD/ AMAZE, but also because of the classification 
scheme, which mainl y relies on a kinemetry analys i s and on 
a visual inspect ion ( |Forster Schreiber et a l. 2009). Forster 
Schreib er et al.| were able to unam biguous ly id entify 18 ro- 
tating systems (discussed in |Cresci et al.|2 009), i.e. a lower 
limit of 35% (18 out of 52 detected systems). These au- 
thors were also able to identify 14 dispersion-dominated 
systems (27% of the SINS sample) based on V max /cr ra- 
tio. Considering only the reliably classified systems, SINS 
and MASSIV percentages are similar. At lower redshift, 
the results of the IMAGES sample (z ~ 0.6) show that 
63% of the sample show signs of rot ation (both rotatin g 
disks and perturbed rotators classes, |Puech et al. 2008). 
This percentage is higher than for MASSIV and favors an 
interpretation in which gas in star-forming systems is stabi- 
lizing into disks while the Universe evolves. This conclusion 
has to be balanced with the various selection functions and 
classification methods, in the same way as for the compari- 
son with the other surveys discussed previously. Concerning 
the selection, the various authors claim that they observed 
representative sets of star-forming galaxies in each redshift 
range. 

In the MASSIV sample, the most massive galaxies 
(logM star > 10.5) are mainly identified as disks in rota- 
tion. However, this is true for the gaseous phase and does 
not necessarily imply that the stars are settled into a disk. 
The best spheroid candidates may be the most compact 
and roundest massive galaxies. Indeed, elliptical galaxies 
can be flattened, too, but should be rounder than disks on 
average. There are four potential candidates even if only 
one of them has an axis ratio fully compatible with zero 
(see Figure 4|. These may also be nearly face-on disks if 
observed velocity shear is very low. In addition, given our 
definition of rotators, it is highly probable that gas and 
stars share a common disk since both morphological and 
kinematic major axis agree well. The morphology obtained 
from the CFHT imaging used for MASSIV has an insuffi- 
cient low spatial resolution to allow addressing this prob- 
lem unambiguously but it is worth noticing that in the local 
Universe stellar disks ar e observed in elliptical galaxies (e.g. 
Kuntschner et al.||2010 ). Usually they are quite young (~ 1 
Gyr) but they can also be as old as ~ 10 Gyr, which would 
be compatible with our observations. 

The intermediate mass galaxies (log M star < 10.5) are 
almost equally divided between rotating and non-rotating 
systems but the smallest ones are principally non-rotating 
systems (see Figure Since both rotating and non- 
rotating systems are observed with small masses and radii, 
we can tentatively conclude that non-rotating systems are 
not caused by observational artifacts. However, we cannot 
exclude the possibility that for some of the latter systems, 
when Rust/ Seeing ~ 0.5, the absence of observed rotation 
is due to an episode of star formation in single non-resolved 
regions. The fact that we observe large non-rotating galax- 
ies incompatible with face-on systems is also very intrigu- 
ing. In addition, those with R e > 2.8 kpc all interact. This 
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Fig. 4. Top: Axis ratio as a function of the stellar mass. 
Bottom: Half-light radius as a function of the stellar mass. 
Same symbols as Figure [3] Arrows indicate that uncertain- 
ties are exceeding the displayed range. 



property is probably related to the nature itself of these 
systems. 

5.2. Dynamical support of disks 

To quantify the dynamical support, it is common to study 
the ratio of the rotation velocity over the local velocity dis- 
persion. In Figure [5j we show this ratio as a function of the 
half-light radius. Only two out of the 20 secure rotators in 
MASSIV have a ratio lower than unity and are consequently 
dispersion-dominated. Only five out of these 20 secure rota- 
tors have V max /(j < 2, which indicates that the majority of 
rotating disks are clearly rotation-supported at z < 1.5. 
Moreover, there seems to exist a trend that the largest 
galaxies have the highest V max /a ratio. As seen in Table |5| 
the median velocity dispersion is a = 62 km s _1 for rota- 
tors. In addition, the median rotation velocity V ma x = 141 
km s _1 and the median V max /<j ratio is 3.2. By restricting 
the analysis on rotating systems with R e > 3.2 kpc (me- 
dian for all rotators), we end up with ten objects with a 
median velocity dispersion a = 49 km s _1 , a median rota- 
tion velocity V = 201 km s _1 and a median V max /a = 3.7, 
which clearly indicates that the largest disks are more sta- 
ble. This seems to be supported in Figure [6] where we clearly 
see that the velocity dispersion is low for the largest rota- 
tors (R e > 6 kpc). 




4 6 

Re [^C] 

Fig. 5. V max /<j as a function of the half light radius. Same 
symbols as Figure [4] 



At higher redshift, rotators in the LSD/AMAZE sam- 
ple (z ~ 3.3) have a typical gaseou s velocity dispersion o f 
a = 90 km s _1 and V max /(j ~ 1.6 ( Gnerucci et al.[ 2011 ), 
whereas in the SINS sample (z ~ 2.2), |Cresci et al.f ( [2009 ) 



found a = 52 km s _1 and V max /a ~ 4.5 as median val 
ues for their unambiguous rotators. The comparison with 
SINS is not straightforward since they have made these 
measurements only for a fraction of these systems (12 out 
of 18). For their sample of rotating systems, the median 
rotation velocity is V max ~ 240 km s _1 and the median 
half-light radius (using the conversion between half-light 
radius and disk scale length R e = l.QSRd valid for an ex- 
ponential disk distribution) is R e = 3.4 kpc. These median 
values are higher than those we obtain for our sample of 
20 rotating galaxies (Vmax ~ 140 km s _1 and R e = 3.2 
kpc), which indicates that while their sample of rotation- 
dominated galaxies probes the massive disks population at 
z ~ 2.2, our sample of rotating systems spans a wider range 
in terms of dynamical mass. At lower redshift, from the 
IMAGES sample (z ~ 0.6), the typical gaseous velocity 
dispersion is a ~ 45 km s _1 (Puech, private communica- 
tion) and Vmax/a ~ 4 considering both rotating disks and 
perturbed ro tators, with a tren d for higher V m ax/ (J f° r ro- 
tating disks ( |Puech et al.|2007[ ). The GHASP sample is the 
largest 2D kinematics sample of local spiral galaxies with 
star-formation rates typical of the local Universe. We con- 
sider here the subsample of 136 galaxies that we projected 
at z ~ 1.33 with similar spati al reso lution conditions as the 
MASSIV sample (see section 3.2.3). These galaxies have a 
median optical radius of 8.5 kpc. Since R opt ~ 1.9R e for an 
exponential disk distribution, this gives a median half-light 
radius of 4.5 kpc. The median rotational velocity is 164 
km s _1 , the median velocity dispersion is 24 km s _1 and 
the median Vmax I & is ~ 7. Despite the small differences 
observed from z ~ 2.2 to z ~ 0.6, observations of various 
samples suggest a trend of decreasing velocity dispersion 
with decreasing redshift. Associated to local and z ~ 3 ob- 
servations, this supports the idea that the gaseous phase of 
galaxy disks becomes less turbulent with cosmic time. 

We also observe that about half of the sample of rotators 
have a gaseous velocity dispersion unambiguously higher 
than 60 km s _1 . These may be interpreted as clumpy disks. 
Indeed, it has been claimed that a high- velocity dispersion 
is expected for clumpy disks created from a smooth cold gas 
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Fig. 6. Velocity dispersion as a function of the effective 
radius. Same symbols as Figure [4] 
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Fig. 7. Velocity dispersion as a function of the redshift. 
Same symbols as Figure [4j The dotted line represents the 
velocity dispersion corresponding to the spectral PSF of 
SINFONI. 



accretion (e.g. Bournaud & Elmegreen 12009] ). This m ode of 
accretion i s believed to be more efficient at z > 2 ( |Keres| 
et al.|2009 ) and it is therefore not surprising to also observe 
a significant percentage of galaxies in MASSIV that resem- 
ble stable rotating disks with velocity dispersions compat- 
ible with - 20 - 50 km s _1 (cf. figure 0. 

We also point out that the gaseous velocity dispersion 
is not significantly different when considering non-rotating 
systems, which raises the question of the origin of this 
gaseous velocity dispersion. If these non-rotating systems 
are indeed transient mergers of star-forming disks in an 
unstable phase, then a high- velocity dispersion is expected 
( Bournaud et al.||2011 ), but then it is more difficult to ex- 
plain cases of large non-rotating galaxies with relatively low 
gaseous velocity dispersion. 

As poin ted out by e.g. Lehnert et al. ( 2009 ) , Green et ah] 
(2010), and Le Tiran et aE" p011| ), star formation might be 
responsible for the gaseous turbulence. In particular, star 
formation intensity might be correlated with gaseous ve- 
locity dispersion. Within this hypothesis, the lowering of 
the velocity dispersion from high to low redshift might be 
explained by the fact that the samples discussed in this pa- 
per have on average a decreasing star formation rate from 



z ~ 3 to z = (see Contini et al. 2012). This behavior is 
due to cosmological surface brightness dimming but also to 
galaxy evolution itself: on average, galaxies have a larger 
size and form less stars at z = than at z ~ 2. The rela- 
tions between star formation and velocity dispersion will be 
addressed with the full MASSIV sample in a forthcoming 
paper. 

5.3. Importance of merging at z ^ 1.2 

The proportion of interacting galaxies in our sample is at 
least 29% (13 interacting systems with flag "A" or "B" out 
of the 45 detected systems with 0.9 < z < 1.6) and is 32% if 
we only consider the 41 systems with a flag "A" or "B" in the 
redshift range 0.9 < z < 1.6. Because of the way we deter- 
mined which systems are in interaction, this proportion is a 
lower limit. Indeed, on the one hand, it could be that a frac- 
tion of non-rotating galaxies are ongoing mergers or merger 
remnants. On the other hand, close mergers could mimic 
rotators (elongation and velocity gradient along a common 
axis), as is the case for VVDS020294045 (see Appendix [A}. 
In addition, owing to the non-circular final field-of-view of 
our SINFONI data we only cover a fraction of the area at a 
given radius around our targets. There is also a substantial 
percentage of objects for which close galaxies are observed 
in I-band but not in SINFONI datacubes (see Appendix 
|A| ). It is possible that a fraction of these objects are passive 
galaxies at the same redshift as the MASSIV targets. It is 
interesting to notice that on average, galaxies in interaction 
are larger than isolated galaxies and have a lower gaseous 
velocity dispersion (see Table [5|. 

The proportion of interacting systems in MASSIV is 
comparable to the percentage of mergers in the SINS sam- 
ple. However, SINS mergers (around one third of the sam- 
ple) are mainly identified using kinemetry, a technique 
based on the degree of perturbation observed in the kine- 
matic maps, whereas the 29% of interacting galaxies in 
MASSIV are mainly based on the detection of several com- 
ponents. On the one hand, only three pairs out of 52 de- 
tected systems were identified in SINS (and three pairs 
out of 29 systems in the LSD/ AMAZE sample), which is 
much lower than the pair fraction observed in MASSIV. 
On the other hand, SINS mergers found from kinemetry 
may m ainly correspond to our non-rotating galaxies (see 
Section 4.4). That means that using this latter method, we 
would probably find that some MASSIV isolated galaxies 
are likely mergers. This is additionally supported by the 
fact that the percentage of non-rotating systems is larger 
in interacting galaxies: among the 11 systems with inter- 
actions in the sample (with both S/N>5 and isolation flag 
"A" or "B"), seven have at least one non-rotating component 
(64%), whereas among the 22 isolated systems (with both 
S/N>5 and isolation flag "A" or "B"), there are only eight 
non-rotating galaxies (36%). This could indicate that some 
isolated non-rotating galaxies are in fact ongoing mergers 
or merger remnants. This seems to be also supported by 
the number of interacting galaxies, which is larger at z ~ 1 
(e.g. Figure [7]). The conclusion is therefore that the per- 
centage of mergers at z ~ 1.2 is higher than at z ~ 2.2 as 
observed in SINS. 

At lower redshift, 44% of galaxies in the IMAGES sam- 
ple have complex kinematics that could be explained by on - 
going mergers of star- forming galaxies (Yang et al. 2008). 
This is more than the proportion of interacting galaxies 
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in MASSIV. However, the definition of galaxies with com- 
plex kinematics is probably closer to our definition of non- 
rotating galaxies (as is the case for the comparison with 
SINS) and the proportions of these two classes in each sam- 
ple are fully comparable. It is probable that for IMAGES, a 
larger field of view would have led to the detection of close 
companions. Consequently, a direct comparison of merging 
with MASSIV is not straightforward. 

The number of mergers deduced from the pair number 
is higher in MASSIV than in SINS. Our results are still 
consistent with a peak in the merger activity at 1 < z < 2 
if we consider systems with complex kinematics. However, 
it is not yet clear if these kinematic signatures of mergers 
are related to minor or major mergers, whereas from the 
MASSIV sample we have a way to infer the mass ratio 
between pair members. A dedicated analysis on the merger 
rate from the observed pairs will be performed for the entire 
MASSIV sample (Lopez-Sanjuan et al., in prep.) and will 
probably enable us to infer the rate of observed ongoing 
mergers. 



Merger remnants in a transient unstable state for the 
gaseous phase (e.g. collapse of gas due to interaction, 
anti-spin mergers, etc.). 

Rotators with a special gaseous distribution: the most 
plausible being a unique highly concentrated star- 
forming region when galaxies are small. 
Galaxies with an unstable gaseous phase. This could 
be supported by the relatively low mass (thus a low 
potential well) of these systems. 

Nearly face-on disks (for low-velocity shear galaxies 
only) : this would explain the similar velocity dispersion 
as for the rotators but the number of these galaxies is 
statistically low. 

Spheroids: however, (i) these galaxies may be the most 
massive ones, (ii) there is no reason to think that the 
gas is not in a disk as observed at lower red shift (e.g. 
Pele tier et al.|[2QQ7| in early- type spirals; e.g. |Bournaud| 



et al.|2QQ7 in numerical simulations of gas-rich mergers) 



and (hi) if the gas is randomly distributed, a higher 
velocity dispersion than for rotators would be expected. 



5.4. Nature of non-rotating galaxies 

Clearly, non-rotating systems are mainly galaxies classified 
as low-velocity shear galaxies plus some interacting galax- 
ies. Among the 16 (and eight additional with S/N < 5) 
low-velocity shear galaxies only three (and two additional 
that have very poor S/N) could be considered as rotating. 
Two of these galaxies are compatible with nearly face-on 
systems. The exact nature of these non-rotating objects 
is still unclear. Such a population of galaxies has already 
been observed at higher redshift both in the SINS sample 
at z - 2.2 and in the LSD/AMAZE sample at z - 3.3. 
These objects are smaller on average than rotators and are 
often associated to interacting systems (see Table [5]). 

An unexpected trend is observed for these objects in the 
MASSIV sample (Figure [6]): there exists an anti-correlation 
between the mean velocity dispersion of the gas and the ef- 
fective radius of the stellar component. This correlation is 
in contrast to what would be expecte d for elliptical galax- 



ies located in the f undamental plane ( |Dressler et al.||1987 
Djorgovski fc Davis|1987 ). However, the fundamental plane 



applies for the central stellar velocity dispersion, which 
is a quantity that cannot currently be measured in high- 
redshift galaxies from spectroscopic data. Alternatively, 
this correlation could reflect a beam-smearing effect since 
an unresolved velocity gradient could reproduce this trend. 
However, the fact that we also observe velocity gradients 
in other galaxies with similar sizes contradicts this inter- 
pretation but it is closely linked to the spatial distribution 
of the ionized gas. For face-on disks (possible only for low- 
velocity shear galaxies), it would be expected to recover the 
same trend as for rotators (i.e. lower velocity dispersion for 
large disks). This is very unlikely, however, because veloc- 
ity gradients may be observed for the largest galaxies and 
also because we do not expect to observe more than four 
galaxies with an inclination lower than 20° (see Figure [2]). 
If the gas were more concentrated, though, it might col- 
lapse more violently to form stars at a higher density if we 
consider that star formation c an drive the gaseous velocity 
dispersion as suggested by e.g. Green et al. (2010). 



As was been proposed in the previous sections, these 
non-rotating objects can be interpreted as 



6. Summary and conclusions 

We presented the data and their processing for a su bsample 
of 50 galaxi es of the MASSIV sample (described in Contini 
|et al.|[2012[ ). We focused on the analysis of the dynami- 
cal state of galaxies using kinematic maps derived from 
SINFONI data-cubes. We classified this sample based on 
the modeling of these maps and based on a comparison 
between morphology and kinematics. This classification de- 
scribes on the one hand the rotating or non-rotating nature 
of the objects and on the other hand the close environment 
of galaxies. 

We separated rotating galaxies and non-rotating galax- 
ies based on two arguments: (i) the agreement between mor- 
phological and kinematic major axis position angle and (ii) 
the accuracy of the description of a rotating disk model. 
We showed that about half of the sample displays rotation, 
whereas one third does not show a dominant ordered ro- 
tation (the remaining fraction is not classified). However, 
the ionized gas turbulence in these two classes is similar 
(a ~ 60 km s _1 ), which marks a transition between higher 
redshift galaxies showing on average a higher velocity dis- 
persion ((j no 60 — 90 km s _1 ), as in SINS at z no 2.2 or 
LSD/AMAZE at z ~ 3.3, and lower galaxies in IMAGES 
at z no 0.6 and GHASP at z = 0, which are characterized 
by a lower velocity dispersion (<j ~ 20 — 40 km s _1 ). All 
these samples are 3D spectroscopic samples observing ion- 
ized gas emission lines. This transition of the gaseous veloc- 
ity dispersion seen from these high- to low-redshift samples 
is partly responsible for the increase of the rotational sup- 
port (deduced from the Vm ax /(j ratio) when galaxies evolve. 
It could be related to a common process that would induce 
a decrease of the star formation rate. We observe that about 
half of our galaxies have a velocity dispersion compatible 
with that in local star-forming galaxies, whereas the other 
half is more compatible with galaxies at higher redshift. At 
high redshift, this high- velocity dispersion seen in isolated 
disks is thought to be caused by a mass assembly driven 
by cold flows. In this framework, our results suggest that 
at z no 1.2 cold gas accretion is less efficient than at higher 
redshift but more efficient tha n at lower redshift, consistent 
with cosmological simulations ( Keres et al.|20 05). When we 
compare our sample to these other high-redshift samples, 



14 



B. Epinat et al.: MASSIV: Mass Assembly Survey with SINFONI in VVDS 



it seems that the fraction of disks increases in star-forming 
galaxies while the Universe evolves. 

By studying strong kinematic signatures of merging and 
detecting pairs in our data-cubes and broad band images, 
we have shown that the fraction of interacting galaxies is 
up to at least one third of the sample. The fraction of non- 
rotating objects in these systems is higher than in isolated 
ones. This suggests that a significant fraction of isolated 
non-rotating objects could be mergers in a transient state 
in which the gas is not dynamically stable. However, the 
nature of non-rotators is still unclear. Compared to higher 
and lower redshift 3D spectroscopic surveys, our findings 
seem to indicate that at the peak of star formation activity, 
the fraction of star-forming galaxies in interaction is also 
at its m aximum, corresponding to a peak in the merging 
activity ([Ryan et al. 2008 Conselice et al. 2008 Lopez 
Sanjuan et al.||2(DQ9 



All these results favor a scenario in which the mass 
assembly of star-forming galaxies is progressively shifting 
from a predominance of smooth cold gas accretion to a 
predominance of merging as cosmic time evolves, with a 
transition epoch around a redshift z ~ 1 — 1.5. 
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Appendix A: Individual comments 

VVDS020106882 This galaxy has a regular velocity 
field. The velocity dispersion is higher on the southeast 
side but remains low. Some external pixels are detected 
on the northwest side. They seem to be related to regions 
associated to the galaxy due to concordant velocities. Two 
objects (including a very faint one) are detected in I-band 
at less than 6" but are not detected in the SINFONI 
datacube. 
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VVDS0201 16027 This galaxy has a small compan- 
ion at 4" to the north. The companion is detected in 
Ha but with a low S/N and over a region of the size 
of the seeing. This companion is less exposed than the 
main galaxy. The main galaxy itself seems to have two 
components in the Ha flux distribution and has a very low 
velocity shear. 

VVDS020147106 Despite a low-velocity gradient 
and some perturbations in the low S/N regions, the 
kinematic position angle of this galaxy agrees well with 
the I-band morphology. This galaxy may therefore be a 
nearly face-on disk (at least the ionized gas) but with a 
high velocity dispersion. There are several objects detected 
in I-band but not in the SINFONI datacube at less than 6". 

VVDS020149061 The I-band image suggests small 
and extremely faint companions that are not detected in 
the SINFONI datacube, however . The seeing is half the 
size of the Ha emission. The velocity field is irregular 
on the edges where the S/N is lower and the velocity 
dispersion is higher than 70 km s _1 everywhere in the 
galaxy (except at the edges). 

VVDS020164388 Except for the edge, the velocity 
field and the velocity dispersion map are regular. Velocity 
dispersion is not higher than 60 km s _1 . In the I-band 
image, the galaxy is quite round. Several faint close objects 
are also observed in this image but are not detected in Ha. 

VVDS020167131 Some [Om] is detected in the center 
but there are only a few pixels and the S/N is lower than 
4. In addition, some emission line is detected at the same 
redshift in the center of the galaxy on the east seen in 
the I-band. This indicates that these two galaxy are in 
interaction. 

VVDS020182331 The velocity gradient is clear even if 
we exclude low S/N regions (bluest and reddest velocities). 
However, the velocity field and velocity dispersion map 
look perturbed, probably because of the low S/N of the 
observation. A very close companion is detected in the 
I-band image that is undetected in Ha. 

VVDS020193070 Despite a generally low S/N, the 
velocity field is quite regular except in the outskirts. The 
velocity dispersion is low but slightly asymmetric. The 
morphology is elongated both in I-band and in Ha (but 
with a difference of ~ 10° between morphological and 
kinematic major axis position angles). Three very faint 
objects are detected at less than 5" from the galaxy but 
are not detected in Ha. 

VVDS020208482 The S/N is fairly low but there 
is a clear velocity shear. The I-band image shows a close 
galaxy and other objects, farther away, that are not 
detected in Ha. 

VVDS020214655 This galaxy displays a fairly smooth 
low velocity gradient, but not aligned with the morphology, 
which is quite round. There are some perturbations in 
the low S/N regions of the velocity field. There are also 
small perturbations in the velocity dispersion map, which 



is fairly flat and has low values. There is a close galaxy in 
the I-band (at ~ 2") that is not detected in Ha. 

VVDS020239133 The morphology from I-band is 
clearly elongated and is more extended than that in Ha. 
The velocity field is asymmetric but shows a clear gradient. 
The velocity dispersion map is also perturbed and has 
a high-velocity dispersion in the south that could be 
explained by a sky line residual. Two objects are observed 
in I-band but not detected in Ha within the SINFONI 
field-of-view. 

VVDS020240675 The velocity field is perturbed, 
with a low- velocity gradient. The velocity dispersion is 
fairly low (~ 40 km s _1 ). In addition, the size of the Ha 
emission is not large compared to the seeing. There is some 
emission detected at 2" in the east of the main galaxy, but 
nothing is detected in the I-band at this location. There 
are also several objects in I-band within 6" that are not 
detected in Ha. 

VVDS020255799 The signal in Ha for this galaxy 
is not extended and with a mean S/N lower than 5. There 
is no clear regular velocity gradient. One object is detected 
at 5" in I-band but not in Ha. 

VVDS020261328 The velocity field is not completely 
regular in the blue and red sides. There is also a peak in 
the velocity dispersion where the S/N is lower. Several 
objects are detected in I-band but not in Ha. 

VVDS020278667 Some Ha is detected, but the ex- 
tent is lower than the seeing, and the S/N is lower than 
5. Therefore its classification is very doubtful. There is a 
close galaxy (~ 1.5") in the I-band but it is not detected 
in Ha. 

VVDS020283083 This galaxy looks asymmetric in 
the I-band and this coincides with the Ha distribution 
and with the main perturbation of the velocity field. This 
could be related to a tidal tail. That is why this galaxy is 
classified as non-isolated. 

VVDS020283830 In the Ha flux map there is a 
faint detection smaller than the seeing that coincides with 
a detection in the I-band (North-East). For the main 
object, the distribution is asymmetric. This is also true in 
the Ha distribution, which is clumpy. These clumps may 
be either star-forming regions in one single rotating galaxy 
or two small objects in the process of merging. 

VVDS020294045 The total velocity shear is larger 
than 50 km s -1 . However, both the I-band and the Ha 
images are consistent with a system involving two compo- 
nents: one main component with a low-velocity gradient 
(south) and one small component with a size comparable 
to the seeing (north). This interpretation is supported by 
the high-velocity dispersion between the two components 
due to beam smearing between the two components at 
different radial velocities. There is a sky line residual at 
the velocity of the small component, but the flux detected 
is noticeable. It is therefore not possible to classify the 
small component. This system is classified as interacting 
with a flag B since it could also be a single object but with 
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strong perturbations. 

VVDS020363717 The I-band shows a close (~ 1"), 
very diffuse and faint object that is not detected in the Ha 
(southwest). The velocity shear is low. The velocity field is 
perturbed in the low S/N regions. The velocity dispersion 
map is also perturbed and the velocity dispersion is high 
(^ 90 km s _1 ) after beam smearing correction. 

VVDS020370467 Both velocity field and velocity 
dispersion map are perturbed. The velocity gradient of 
the model is low and the velocity dispersion is quite high 
(higher than 70 km s _1 ) in particular at the north edge 
(higher than 150 km s _1 ). This could be a sign for a strong 
interaction due to merging, but there is also a sky line 
residual that could induce the line width. In addition some 
objects are observed in the I-band image at less than 6" 
but are not detected in Ha. 

VVDS020386743 The Ha distribution is very elon- 
gated for this galaxy, but no velocity gradient is visible. 
Both velocity field and velocity dispersion are smooth. The 
I-band image shows a distorted galaxy compatible with the 
Ha distribution. In the north, the distortion is compatible 
with a faint detection in Ha that could be related to a 
minor close companion. 

VVDS020461235 The velocity field of this galaxy 
is perturbed. There is a blob detected in Ha in the 
southwest side. This blob is also suggested in the I-band 
morphology but is less clear. This blob has velocities 
compatible with the rotation of the main component but 
could as well be a minor companion in the process of 
merging. This system is therefore classified as interacting 
with a flag B. 

VVDS020461893 The velocity field is regular but 
the velocity dispersion is less regular and high and could 
suggest a high-dispersion clump. An object is observed at 
3" southward in the I-band image but is not detected in 
SINFONI data. 

VVDS020465775 Both velocity field and velocity 
dispersion map are irregular. The gradient of the model is 
lower than 50 km s _1 . We claim that there is a companion 
in the northwest that induces broad and non- Gaussian 
(hence a lower S/N) lines interpreted as a sign of interac- 
tion that could also be responsible for the asymmetry in 
the I-band. Hence it is classified as interacting with a flag 
B since it could also be a single object but with strong 
perturbations. 

VVDS140083410 The velocity field is perturbed 
and there is no clear overall velocity gradient. No galaxy is 
detected at less than 8" in the I-band image. 

VVDS140096645 The velocity field is smooth but 
with some perturbations along the minor kinematic axis. 
The velocity dispersion map is perturbed and larger 
on the southern side probably due partially to a sky 
line residual. Some emission is detected in the Ha map 
outside the main component (north). This detection is 
smaller than the seeing and is not clearly confirmed in the 
I-band morphology (CFHT12k images are less deep than 



CFHTLS images). It could be a minor companion but it is 
also compatible with the velocity of the main component 
and could be associated to it. Therefore it is classified as 
intera cting with a flag B. This galaxy also has an AGN 
see IQueyrel et aL|2012| >. 

VVDS140123568 Some Ha is detected, but there 
are only a few pixels and the S/N is lower than 4. 
Therefore its classification is very doubtful. There is 
no detected object in the I-band image at less than 6" 
around the galaxy. Since the kinematics does not allow to 
investigate a possible ongoing merger, the isolation has a 
flag B. 

VVDS140137235 Some Ha is detected, but there 
are only a few pixels and the S/N is lower than 4. 
Therefore its classification is very doubtful. There is 
no detected object in the I-band image at less than 6" 
around the galaxy. Since the kinematics does not allow to 
investigate a possible ongoing merger, the isolation has a 
flag B. 

VVDS140217425 The velocity field is regular and 
has the highest shear (~ 600 km s _1 ). The velocity disper- 
sion map shows two regions with high-velocity dispersions 
(> 150 km s -1 ). These regions are associated with a lower 
S/N. This is because there are double profiles in these 
regions that induce a poor fit. These double profiles are 
caused by the beam-smearing that mixes regions with 
different velocities. This can be interpreted as the presence 
of star-forming clumps at the edge of the galaxy and a 
large clump in the center. The decomposition of the profile 
using two Gaussian scales down the velocity dispersion to 
usual values, and shows that the two external clumps have 
no velocity gradient, whereas the central one does have 
one. This picture is compatible with a rotating disk in 
which the rotational velocity reaches a plateau before the 
external clumps. Using the double profile decomposition 
instead of the kinematics model to derive the maximum 
velocity, we find V max ~ 312± ~ 20 km s _1 , which is much 
more realistic than the value found from the model. The 
velocities of the plateau are ~ +313 and ~ —312, which 
shows that the rotation curve is fairly symmetric and 
therefore favors a rotating disk hypothesis. Correcting for 
the inclination, we obtain V max = 322 km s _1 . However, 
the I-band image is asymmetric (which is not incompatible 
with the previous statement), and it cannot be ruled out 
that the external clumps are instead objects about to 
merge with the main component. 

VVDS140258511 The velocity field and velocity 
dispersion map look perturbed at the northeast edge 
but this may be due to a sky line residual. Except at 
this location, the velocity dispersion is low (less than 40 
km s _1 ). 

VVDS140262766 The velocity field and the velocity 
dispersion map are not strongly perturbed and the velocity 
dispersion is about 50 km s _1 . Several objects are seen at 
less than 5" in the I-band image but are not detected in Ha. 

VVDS140545062 This object has a clear and smooth 
velocity shear compatible with the I-band and Ha mor- 
phology. However, the velocity field is not very well 



B. Epinat et al.: MASSIV: Mass Assembly Survey with SINFONI in VVDS 



17 



reproduced by the rotating disk model: one side is faster 
than the other. The velocity dispersion map shows a peak 
in the southwest side. This peak is caused by a double 
profile that may be interpreted as a sign of merging (late 
stage or ongoing merger). However, since the system as a 
whole is isolated (no galaxy at less than 1" in I-band), this 
galaxy was classified as isolated with a flag B. 

VVDS220014252 Both velocity field and velocity 
dispersion map are perturbed. In addition, the morphology 
in I-band is asymmetric. This system has a high- velocity 
dispersion (> 80 km s _1 ) and is peaked on some edges. 
This may be a sign of merging (ongoing or remnant). 
However, the system is isolated and there is no other evi- 
dence for some companion except one very faint object in 
I-band at ~ 5" which is not in the SINFONI field- of- view. 
So this system is classified as isolated. 

VVDS220015726 The velocity field of this galaxy 
is well reproduced by a rotating disk model. The velocity 
dispersion map is peaked in the center, which agrees with 
the expected effect of beam smearing caused by the inner 
velocity gradient. The velocity field and velocity dispersion 
map are slightly perturbed. One very faint component 
is detected in I-band at around 4" but is not detected in Ha. 

VVDS220148046 This galaxy has been observed 
using AO and is therefore less deep than seeing-limited 
observations. Some emission line was detected but with a 
low S/N (~ 4). It appears that the redshift determined 
from the VVDS was wrong and that we observed [Om] at 
z = 2.2442. The emission is compact and therefore it is not 
possible to conclude about the kinematics. 

VVDS220376206 A small component (smaller than 
the seeing, hence a flag B for the environment classifica- 
tion) is detected in Ha and is confirmed in the I-band 
image (north). The residual velocity field shows signs of 
interaction. The main component has two blobs in Ha 
and the I-band morphology is slightly asymmetric. The 
velocity dispersion map is also perturbed and peaked on 
the edges. This may be a sign of interaction or of the end 
of a merging event. 

VVDS220386469 This observation was made with 
AO and is underexposed (owing to the small pixel size, 
cf. |Contini et al. |2Q12[ ). Even if the AO observation 



does not allow a fair comparison with other galaxies 
in the same conditions because of the smaller field of 
view, this system is classified as isolated with a flag B 
since the closest clear detection in the I-band image is 
farther than the 6" that would be reached using the 
seeing-limited observing strategy. The velocity field is 
irregular and the velocity dispersion is very small but 
this may be due to the use of AO, which provides more 
details and a lower S/N than seeing-limited observations. 
The extent of the velocity field is small so it is difficult to 
argue that there is no velocity gradient at all in this galaxy. 

VVDS220397579 This system is the clearest inter- 
acting system between two large galaxies. The two objects 
are detected in Ha and in the I-band. Both are elongated 
toward the same direction. The main object in Ha and 
in I-band has a smooth velocity field with no gradient 



except on the northern side, near the companion. This 
region coincides with a higher velocity dispersion and this 
is probably related to the interaction with the compan- 
ion. Since the galaxies are aligned and since the overall 
velocity gradient is monotonic, the whole system might 
be considered as a unique galaxy with two big clumps. 
However, such a system would be very much extended, a 
property which is physically unlikely at these redshift s. 
The companion shows a clear velocity shear and has 
been classified as well. It is also worth noticing that 
the minor companion in I-band is the most luminous in 
the K-band (from UKIDSS Deep Extragalactic Survey, 
Lawrence et al. 2007), therefore likely the most massive one. 



VVDS220544103 This galaxy has a clear velocity 
gradient. The object is distorted in the I-band image as 
well as in Ha. This can be explained by the presence of two 
components (the morphology was decomposed using this 
hypothesis). In that case, the southern component may 
be well described by a rotating galaxy and the brightest 
blob in Ha (north) maybe a companion in a merging stage 
or a relic of merging. However, it cannot be ruled out 
that this is a single object with disturbed morphology and 
kinematics, hence the flag B. 

VVDS220544394 The Ha map clearly shows the 
presence of a main component and a faintest one in the 
north. Their velocity fields seem decorrelated enough to 
claim that these are two separated objects on a pre-merging 
stage. The small component is also clearly detected in 
the I-band. The velocity field of the main component 
is perturbed but shows a velocity shear. Its velocity 
dispersion is about 50 km s _1 . The small component is 
slightly larger than the seeing and may also be considered 
to have a velocity shear. 

VVDS220576226 The velocity field is perturbed in 
the north where the S/N is the lowest. The velocity 
dispersion map is regular. Two objects are observed within 
6" in the I-band image but are not detected in Ha. 

VVDS220578040 The I-band image shows an asymmet- 
ric galaxy with a more diffuse emission in the east. The Ha 
map shows two blobs. These can be interpreted either as 
two distinct galaxies or as two clumps in one single galaxy. 
The velocity field is perturbed and is not fully reproduced 
by the model. The velocity dispersion map is slightly 
perturbed. Since no strong kinematics perturbation is 
observed, this galaxy is classified as isolated (with a flag 
C) and therefore its kinematics is interpreted as a sign of 
rotation. 

VVDS220584167 The I-band morphology is distorted 
and the Ha distribution is asymmetric. The kinematic 
position angle, however, agrees very well with the mor- 
phology. The velocity field shows some perturbations 
such as a slight asymmetry. This could be due to the 
position of the kinematic center that would differ from 
the morphological one. The velocity dispersion map is 
peaked in the center, however, but the beam smearing 
may not account for the whole amplitude of the peak. Two 
objects are detected in I-band but not in Ha at less than 6". 

VVDS220596913 This galaxy is a chain galaxy in 
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the I-ba nd. It has been obse rved both with and without 
AO (cf. IContini et all 120121). Several blobs are detected 
in Ha. They can be interpreted as several clumps in a 
single edge-on galaxy since the velocity field is monotonic. 
They can also be interpreted as several small galaxies in 
a merging stage. The continuum was extracted from the 
SINFONI data around Ka and seems to be peaked in the 
center, between the Ha blobs, which favors the edge-on 
hypothesis. Since no strong kinematics perturbation is 
observed, and since a previous seeing-limited observation 



Epinat et al. 2009) revealed that no extra emission is 
etected, this galaxy is classified as isolated with a flag 

B and therefore its kinematics is interpreted as sign of 

rotation. 



VVDS910193711 This observation was made with 
AO. The extent is larger than 1" . The velocity field is 
perturbed. The interpretation of the perturbations may not 
be comparable with seeing-limited observations because 
the S/N is lower and the scale of the perturbations smaller. 
The velocity dispersion is high (around 100 km s _1 ). In 
the close neighborhood that would be observed using the 
seeing-limited strategy, several objects are detected in the 
I-band image. 

VVDS910279515 This galaxy was observed using 
AO and is therefore less deep than seeing-limited observa- 
tions. This galaxy is unclassified since there are only a few 
pixels with a S/N larger than 3. There are several objects 
in the I-band image at less than 5" that would have been 
observed using a seeing-limited strategy. 
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02:26:28.926 


A A A A O A A A 

—04:22:31.14 


1.3949 


8.43 


A -1 A f A A f 

0.125x0.25 


H 


Ha 


120 


0.77 


2472 


P82G 


020294045 


02:25:47.152 


—04:21:07.41 


1.0028 


8.01 


0.125x0.25 


J 


Ha 


120 


0.59 


2067 


P79B 


020306817 1 


02:25:50.316 


-04:19:22.93 


1.2225 


8.32 


0.125x0.25 


J 


[Om] 


120 


— 




P79B 


020363717 


02:26:23.709 


-04:11:57.87 


1.3339 


8.40 


0.125x0.25 


H 


Ha 


80 


0.64 


2 107 


P82G 


020370467 


02:26:14.690 


-04:11:05.44 


1.3338 


8.40 


0.125x0.25 


H 


Ha 


80 


0.71 


2407 


P82G 


020386743 


02:27:13.989 


-04:08:59.73 


1.0487 


8.09 


0.125x0.25 


J 


Ha 


120 


0.73 


2111 


P79B 


020461235 


02:26:47.102 


-04:23:55.70 


1.0349 


8.07 


0.125x0.25 


J 


Ha 


120 


0.63 


2101 


P79B 


020461893 


02:27:12.252 


-04:23:11.28 


1.0486 


8.09 


0.125x0.25 


J 


Ha 


80 


0.60 


2115 


P82G 


020465775 


02:26:59.366 


-04:19:00.08 


1.3583 


8.41 


0.125x0.25 


H 


Ha 


80 


0.88 


2434 


P82G 


140083410 


13:57:50.595 


+04:17:38.71 


0.9435 


7.89 


0.125x0.25 


J 


Ha 


80 


0.69 


2005 


P81D 


140096645 


13:58:26.336 


+04:19:47.75 


0.9655 


7.94 


0.125x0.25 


J 


Ha 


120 


0.56 


2024 


P81D 


140123568 


13:55:57.628 


+04:24:20.11 


1.0012 


8.01 


0.125x0.25 


J 


Ha 


120 


0.76 


2067 


P79A 


-i a r\ a o ^ a o r^ - 

140137235 


-i o r~ r* a a ^ a a 

13:56:12.729 


i A A a r* o -i T A 

+04:26:31.74 


1.0445 


8.09 


A -1 A r* A A P" 

0.125x0.25 


J 


Ha 


80 


0.76 


2111 


P79A 


i a r\c\A <—7 a c\f 

140217425 


13:57:56.405 


+04:38:37.00 


0.9792 


7.97 


0.125x0.25 


J 


Ha 


100 


0.95 


2040 


P81D 


1 A AO C ori 1 

i4Uzoooll 


1 A . AA . 1 A CCO 

14:UU:19.ooo 


i A A . A A . A C OC 

+U4:44:4o.oo 


1 O A OO 

l.z4zo 


o o o 

8.33 


A IOCwA OC 

U.lzoxU.zo 


TT 
11 


T T „ 

Ha 


OA 

8U 


A A A 

U.49 


OO 1 c 

2315 


TDTA A 

P79A 


14Uzoz7oo 


10 C A C C C 1 O 

13:59:55.518 


lA/l AC OA A/I 

+04:45:30.04 


1.2836 


8.37 


A IOCwA OC 

O.lzoxO.zo 


TT 

li 


T T „ 

Ha 


120 


0.51 


2353 


P79A 


1 A AC a cnco 

i4Uo45Uoz 


lo:59:oo.59o 


l A C . O A . O 1 11 

+05:30:31. 11 


1 A A AO 

1.0408 


O AO 

8.08 


A IOCwA OC 

U.lzoxU.zo 


T 
J 


Ha 


1 OA 

IzU 


A TA 

U.7U 


O 1 A A 
Z1U4 


"DO 1 TA 


OOAA1 A O C O 

zzUU14zoz 


OO IT An OTT 


i AA OO OA CO 

+00:28:39. 5z 


1.3105 


8.38 


A 10C\xA OC 

U.lzoxU.zo 


TT 

li 


T T „ 

Ha 


120 


0.70 


2384 


P75A 


OOAA1 C TO/? 

ZZUUlOYZO 


OO. 1 C . /I O /IOC 

zz:io:4z.4oo 


l AA.OA.AO CO 

+00:29:03.58 


1 OAOO 

1.2933 


O OT 

8.37 


A lOCv^A OC 

U.lzoxU.zo 


TT 

li 


Ha 


1 OA 

IzU 


a a a 
U.46 


OOO A 

2384 


T) TATD 

P79B 


220071601 1 


22:18:01.569 


+00:45:34.69 


1.3538 


8.41 


0.050x0.10 


H 


Ha 


80 






P79B 


A A A i /I A A /I /"» 

220148046 


22:14:37.904 


i a -i aa aa r* r~ 

+01:08:20.65 


A A A A A 

2.2442 


8.24 


A A P" A A "1 A 

0.050x0.10 


H 


[Om] 


80 


0.27 


2450 


P81E 


220376206 


22:20:05.772 


—00:08:21.74 


1.2445 


8.34 


0.125x0.25 


H 


Ha 


120 


0.50 


2315 


P79B 


220386469 


22:19:56.603 


-00:03:03.78 


1.0226 


8.05 


0.050x0.10 


J 


Ha 


40 


0.23 2 


2090 


P79B 


OOAOATCTA 


OO.OA.OC CIO 

zz:zU:oo.olz 


1 AA.A1 .Ad OK 

+UU:Ul:4o.oo 


1 A OTA 
1.U3 ( 9 


O AO 

o.Uo 


A IOCwA OC 

U.lzoxU.zo 


J 


TJ 

Ma 


1 OA 

IzU 


a a a 
U.D4 


01 A1 

ZlUl 




220544103 


22:15:25.689 


+00:06:40.31 


1.3973 


8.43 


0.125x0.25 


H 


Ha 


120 


0.76 


2472 


P75A 


220544394 


22:14:24.153 


+00:06:46.67 


1.0101 


8.03 


0.125x0.25 


J 


Ha 


120 


0.58 


2073 


P79B 


220576226 


22:16:11.417 


+00:16:30.46 


1.0217 


8.05 


0.125x0.25 


J 


Ha 


120 


0.58 


2087 


P79B 


220578040 


22:17:04.113 


+00:16:56.80 


1.0462 


8.09 


0.125x0.25 


J 


Ha 


120 


0.62 


2111 


P79B 


220584167 


22:15:22.917 


+00:18:48.82 


1.4655 


8.45 


0.125x0.25 


H 


Ha 


120 


0.75 


2541 


P75A 


220596913 


22:14:29.179 


+00:22:18.93 


1.2658 


8.35 


0.050x0.10 


H 


Ha 


120 


0.18 


2340 


P79B 


910193711 


02:25:46.285 


-04:32:33.43 


1.5564 


8.47 


0.050x0.10 


H 


Ha 


80 


0.27 


2636 


P82F 


910279515 


02:25:36.233 


-04:21:16.13 


1.4013 


8.43 


0.050x0.10 


H 


Ha 


80 


0.21 


2477 


P82F 



The coordinates of the galaxies are given in column (2) and (3). The redshift (4) is determined from SINFONI data. The physical 
scale (5) is computed from the redshift. The SINFONI pixel scale is given in column (6). In column (7), the wavelength band used 
to to observe the emission line (8) is given. Column (9) is the on-source exposure time. The seeing of SINFONI observations (10) 
is derived from observed PSF stars. Column (11) is the spectral resolution. The ESO observing period is given in column (12). 

1 These galaxies are those for which no line was detected in SINFONI data. The redshift is the one derived from VIMOS spectra, 
the line is the targeted one and both spectral and spatial resolution have not been measured. 

2 No PSF star was observed: the resolution is the mean resolution of AO observations. 
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Table 2. CFHT observations and morphological parameters from GALFIT modeling 



VVDS ID Survey Pixel scale Seeing PA b/a R e n 



(1) 


(2) 


["] 

(3) 


["] 

(4) 


[°] 

(5) 






(6) 








[kpc 
(7) 










(8) 






020106882 


LS 





186 





604 


294 ± 4 





62 


± 





03 


3 


51 


± 





14 





77 


zb 





16 


020116027 


LS 





186 





604 


184 ± 2 





38 


± 





03 


4 


27 


± 





16 





70 


zb 





15 


020126402 


LS 





186 





615 


15 ± 8 





56 


± 





10 


2 


12 


zb 





31 





59 


zb 





58 


020147106 


LS 





186 





615 


310 ± 1 





19 


± 





15 


1 


17 


± 





41 





41 


zb 


1 


82 


020149061 


LS 





186 





649 


201 ± 369 





93 


± 





83 


1 


09 


zb 





69 





11 


zb 


6 


68 


020164388 


LS 





186 





618 


162 ± 3 





71 


± 





02 


2 


72 


± 





05 





85 


zb 





09 


020167131 


LS 





186 





618 


272 ± 3 





60 


± 





03 


2 


60 


zb 





09 


1 


64 


zb 





26 


020167131s 


LS 





186 





618 


57 ± 1 





37 


± 





02 


4 


11 


zb 





08 





31 


zb 





09 


020182331 


LS 





186 





618 


268 ± 4 





49 


zb 





04 


4 


19 


zb 





17 





63 


zb 





21 


020193070 


LS 





186 





615 


215 ± 2 





33 


zb 





03 


3 


49 


zb 





09 





47 


zb 





13 


020208482 


LS 





186 





600 


355 ± 7 





70 


± 





05 


3 


67 


zb 





20 





82 


zb 





21 


020214655 


LS 





186 





604 


32 ± 4 





41 


± 





06 


1 


47 


zb 





08 


1 


97 


zb 





53 


020217890 


LS 





186 





604 


171 ± 2 





21 


± 





04 


3 


57 


zb 





14 


1 


13 


zb 





22 


020239133 


LS 





186 





608 


123 ± 2 





39 


± 





03 


3 


14 


zb 





09 





25 


zb 





13 


020240675 


LS 





186 





619 


190 ± 13 





62 


± 





13 


1 


09 


zb 





16 


2 


03 


zb 


1 


19 


020255799 


LS 





186 





607 


164 ± 12 





80 


± 





08 


1 


89 


zb 





26 





24 


zb 





44 


020261328 


LS 





186 





634 


171 ± 8 





51 


± 





11 


1 


83 


zb 





27 





64 


zb 





53 


020278667 


LS 





186 





596 


140 ± 12 





74 


± 





09 


2 


82 


zb 


1 


66 


9 


60 


zb 


5 


68 


020283083 


LS 





186 





596 


301 ± 3 





38 


± 





05 


4 


29 


zb 





20 





46 


zb 





18 


020283830 


LS 





186 





596 


142 ± 1 





35 


± 





01 


6 


82 


zb 





16 





09 


zb 





06 


020294045 


LS 





186 





606 


1 ± 3 





61 


zb 





03 


2 


89 


zb 





10 


1 


00 


zb 





17 


020306817 


LS 





186 





606 


86 ± 2 





54 


± 





02 


4 


22 


zb 





09 





27 


zb 





08 


020363717 


LS 





186 





626 


154 ± 6 





55 


± 





08 





72 


zb 





09 


3 


59 


zb 


1 


01 


020370467 


LS 





186 





626 


50 ± 24 





78 


± 





18 


1 


31 


zb 





25 





54 


zb 





97 


020386743 


LS 





186 





639 


203 ± 8 





64 


± 





07 


2 


68 


zb 





21 


1 


15 


zb 





42 


020461235 


LS 





186 





597 


332 ± 3 





57 


± 





02 


3 


96 


zb 





12 





47 


zb 





11 


020461893 


LS 





186 





634 


283 ± 3 





46 


± 





04 


2 


66 


zb 





11 


1 


04 


zb 





23 


020465775 


LS 





186 





623 


146 ± 3 





51 


± 





03 


4 


04 


zb 





16 





09 


zb 





17 


140083410 


12K 





204 





958 


311 ± 10 





67 


± 





09 


1 


93 


zb 





19 


1 


90 


zb 





82 


140096645 


12K 





204 





747 


217 ± 114 





93 


± 





28 


1 


80 


zb 





40 





85 


zb 


1 


26 


140123568 


12K 





204 





968 


227 ± 229 





80 


± 


1 


14 


1 


39 


zb 





98 


0. 


08 


zb 


32 


.82 


140137235 


12K 





204 





920 


116 ± 24 





47 


± 





31 


5 


45 


zb 


2 


03 





88 


zb 


1 


73 


140217425 


12K 





204 





713 


260 ± 1 





25 


± 





01 


8 


96 


zb 





13 





33 


zb 





04 


140258511 


12K 





204 





768 


265 ± 42 





54 


± 





51 


2 


62 


zb 


2 


11 





40 


zb 


3 


22 


140262766 


12K 





204 





613 


142 ± 28 





65 


± 





26 


2 


07 


zb 





70 





63 


zb 


1 


34 


140545062 


12K 





204 





724 


216 ± 41 





46 


± 





35 


2 


89 


zb 


6 


86 





06 


zb 


4 


03 


220014252 


LS 





186 





687 


136 ± 2 





22 


± 





06 


3 


32 


zb 





12 





94 


zb 





23 


220015726 


LS 





186 





686 


193 ± 12 





78 


± 





07 


2 


76 


zb 





25 





26 


zb 





35 


220071601 


LS 





186 





665 


74 ± 1 





31 


± 





01 


8 


13 


zb 





13 





21 


zb 





04 


220148046 


LS 





186 





800 


236 ± 10 





10 


± 





87 


1 


98 


zb 





38 


2 


10 


zb 


2 


36 


220376206 


LS 





186 





749 


235 ± 1 





46 


± 





01 


5 


29 


zb 





08 





21 


zb 





05 


220386469 


LS 





186 





725 


168 ± 7 





70 


± 





05 


2 


90 


zb 





13 





73 


zb 





26 


220397579 


LS 





186 





707 


331 ± 5 





43 


± 





06 


3 


02 


zb 





32 





21 


zb 





37 


220397579s 


LS 





186 





707 


332 zb 2 





14 


± 





03 


9 


77 


zb 


3 


29 


3 


77 


zb 


1 


49 


220544103 


LS 





186 





686 


208 =b 55 





08 


± 





16 


5. 


61 


zb 


19 


.89 





02 


zb 





91 


220544103s 


LS 





186 





686 


161 ± 5 





26 


± 





06 


5 


47 


zb 





44 


1 


67 


zb 





77 


220544394 


LS 





186 





624 


234 ± 7 





70 


± 





04 


3 


40 


zb 





15 





48 


zb 





13 


220544394s 


LS 





186 





624 


233 =b 10 





30 


± 





17 


3 


69 


zb 





70 


1 


18 


zb 


1 


02 


220576226 


LS 





186 





675 


237 ± 15 





85 


± 





06 


2 


16 


zb 





12 


1 


58 


zb 





36 


220578040 


LS 





186 





658 


112 zb 14 





89 


± 





04 


3 


81 


zb 





12 





42 


zb 





10 


220584167 


LS 





186 





684 


193 ± 2 





52 


zb 





01 


7 


17 


zb 





19 





14 


zb 





06 


220596913 


LS 





186 





623 


255 zb 1 





18 


± 





02 


9. 


49 


zb 


11 


.54 





04 


zb 





12 


910193711 


LS 





186 





597 


9 ± 3 





50 


± 





03 


2 


27 


zb 





06 





50 


zb 





17 


910279515 


LS 





186 





602 


131 ± 3 





61 


± 





02 


4 


44 


zb 





11 





28 


zb 





08 



In column (2), "LS" refers to the CFHT Legacy Surve y (http://www.cfht.hawaii.edu/Science/CFHLS/) whereas "12K" refers to 



the CFH12K/CFHT survey (McCracken et al. 



2003). The pixel scale of the CFHT image and the seeing determined from close 
stars are respectively given in column (3) and (4). The parameters of the model are the position angle of the major axis (5), 
the axis raio (6), the effective radius (7), the Sersic index (7) and the magnitude (not given here due to a zero point magnitude 
mismatch between the two CFHT surveys). 

The suffix "s" refers to secondary objects which have been detected both in Ha and in the I-band image. The companion's 
parameters are only given when they are resolved. 
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Table 3. Physical parameters from kinematic modeling 



VVDS ID 


i 


PA k 




n 




v t 


t r I 

Vt/n 


Resy 


2 

X 


Vmax 




a 




Rlast 




[°] 


[°] 




[kpc] 




n — 1 1 

[km s ] 


n —1 i — li 

[km s kpc J 


n — li 

[km s ] 




n — li 

[km s ] 


[km s 


- 1 ] 


[kpc] 


(1) 


(2) 


(3) 




(4) 




(5) 


(6) 


(?) 


(8) 


(9) 


(10) 


(11) 


020106882 


52 ± 3 


317 zb 1 


1 


1 zb 


1 


133 zb 3 


126 


12 


1.6 


133 zb 25 


41 


zb 


31 


5.1 


020116027 


68 ± 4 


207 zb 6 


4 


3 zb 3 





27 zb 6 


6 


10 


3.7 


27 zb 10 


47 


zb 


15 


6.5 


020147106 


60 zb 24 


317 zb 3 


1 


6 zb 4 


6 


26 zb 3 


16 


6 


1.9 


26 zb 51 


81 


zb 


10 


7.8 


020149061 


60 =b 24 


235 zb 3 




42.0 




976 


23 


15 


1.9 


112 zb 216 


73 


zb 


20 


4.8 


020164388 


45 ± 4 


97 zb 4 


7 


2 zb 4 


9 


79 zb 44 


11 


8 


0.9 


79 zb 19 


52 


zb 


19 


8.2 


020167131 


53 zb 4 


185 zb 10 


1 


zb 


1 


127 zb 43 


123 


13 


0.9 


127 zb 29 


26 


zb 


37 


1.8 


020182331 


61 zb 2 


234 zb 4 




27.0 




647 


24 


11 


0.9 


132 zb 26 


67 


zb 


29 


5.5 


020193070 


71 =b 4 


184 zb 2 




108.7 




3217 


30 


15 


1.8 


117 zb 23 


33 


zb 


27 


3.9 


020208482 


46 zb 3 


336 zb 10 


1 


zb 


1 


158 zb 38 


157 


11 


0.6 


158 zb 31 


7 


zb 


11 


1.4 


020214655 


66 zb 11 


341 zb 4 


1 


zb 


1 


52 zb 4 


51 


8 


1.1 


52 zb 14 


63 


zb 


25 


5.7 


020239133 


67 ± 4 


109 zb 6 




24.6 




766 


31 


13 


0.9 


149 zb 33 


75 


zb 


35 


4.8 


020240675 


60 zb 24 


183 zb 12 




8.9 




116 


13 


13 


1.5 


50 zb 97 


33 


zb 


17 


3.8 


020255799 


37 zb 15 


89 zb 76 


1 


zb 


1 


14 zb 16 


14 


10 


0.6 


14 zb 26 


76 


zb 


23 


4.0 


020261328 


59 ± 11 


179 zb 2 




6.2 




154 


25 


10 


1.6 


127 zb 35 


54 


zb 


19 


5.1 


020278667 


42 zb 20 


175 zb 6 




3.7 




224 


60 


16 


1.5 


77 zb 189 


52 


zb 


37 


1.3 


020283083 


68 zb 2 


359 zb 4 


1 


zb 


1 


59 zb 5 


57 


11 


1.4 


59 zb 12 


39 


zb 


23 


5.6 


020283830 


70 zb 1 


156 zb 2 


1 


9 zb 3 


4 


186 zb 5 


96 


21 


4.7 


186 zb 30 


17 


zb 


24 


7.9 


020294045 


52 zb 5 


3 zb 1 




31.1 




1332 


43 


33 


12.4 


234 zb 51 


60 


zb 


45 


5.5 


020363717 


60 ± 24 


106 zb 6 




16.4 




123 


8 


9 


1.3 


45 zb 86 


91 


zb 


14 


6.0 


020370467 


39 =b 14 


19 zb 11 


1 


zb 


1 


51 zb 10 


49 


18 


1.2 


51 zb 64 


86 


zb 


35 


5.5 


020386743 


50 zb 4 


136 zb 3 




37.0 




286 


8 


9 


2.5 


42 zb 10 


53 


zb 


18 


5.4 


r\c\r\ a r* -i r\ o p" 

020461235 


55 zb 1 


351 zb 2 


1 


2 zb 7 


7 


82 zb 4 


67 


11 


1.8 


82 zb 16 


24 


zb 


22 


5.4 


020461893 


63 zb 4 


279 zb 4 


1 


2 zb 9 


1 


58 zb 6 


46 


8 


1.1 


58 zb 13 


67 


zb 


22 


6.5 


020465775 


59 zb 2 


178 zb 5 


1 


1 zb 


1 


68 zb 7 


65 


13 


1.5 


68 ± 15 


84 


zb 


30 


4.9 


140083410 


48 ± 16 


39 zb 6 


1 


zb 


1 


30 zb 4 


30 


10 


2.3 


30 zb 33 


66 


zb 


21 


5.3 


140096645 


22 zb 15 


197 zb 1 


1 


1 zb 1 


1 


295 zb 12 


258 


12 


6.1 


295 zb 709 


77 


zb 


27 


4.5 


140123568 


60 =b 24 


184 zb 14 




1.0 




108 


108 


6 


0.7 


50 zb 99 


73 


zb 


24 


0.5 


140137235 


62 dz 1 


123 zb 16 




2.2 




200 


91 


4 


0.6 


62 zb 11 


17 


zb 


27 


0.7 


140217425 


76 zb 2 


258 zb 1 




16.1 




499 


31 


29 


20.7 


320 zb 46 1 


45 


zb 


33 


14.5 


140258511 


57 dz 3 


213 zb 1 


1 


zb 


1 


124 zb 4 


119 


19 


5.2 


124 zb 26 


25 


zb 


29 


5.2 


140262766 


60 zb 24 


175 zb 2 




18.1 




524 


29 


9 


1.2 


119 zb 231 


39 


zb 


16 


4.1 


140545062 


63 zb 5 


229 zb 1 




13.9 




378 


27 


13 


3.6 


204 zb 46 


67 


zb 


29 


7.5 


220014252 


77 zb 2 


141 zb 1 


1 


zb 


1 


129 zb 1 


123 


15 


4.5 


129 zb 27 


90 


zb 


28 


10.3 


220015726 


39 zb 15 


186 zb 1 


1 


5 zb 


4 


231 zb 4 


155 


10 


3.1 


231 zb 356 


62 


zb 


22 


3.7 


220148046 


60 zb 24 


261 zb 15 




1.3 




61 


49 


10 


2.9 


42 zb 83 


46 


zb 


21 


0.9 


220376206 


63 zb 2 


225 zb 1 


8 


3 zb 


2 


201 zb 4 


24 


15 


4.8 


201 zb 27 


73 


zb 


26 


10.0 


220386469 


46 zb 7 


151 zb 10 




63.0 




973 


15 


16 


2.0 


40 zb 11 


43 


zb 


25 


2.6 


220397579 


65 zb 6 


1 zb 9 




35.9 




32 


1 


9 


7.0 


9 zb 10 


59 


zb 


17 


10.2 


220397579s 


80 zb 5 


344 zb 2 


3 


8 zb 1 


2 


222 zb 9 


59 


17 


2.6 


222 zb 15 


27 


zb 


37 


6.2 


220544103 


80 zb 4 


198 zb 1 


1 


1 zb 


1 


137 zb 2 


130 


12 


4.4 


137 zb 24 


71 


zb 


19 


7.6 


220544394 


46 zb 2 


180 zb 2 


2 


zb 1 


5 


55 zb 4 


27 


5 


0.8 


55 zb 11 


49 


zb 


17 


5.0 


220544394s 


73 =b 5 


198 zb 4 


1 


zb 


1 


94 =b 8 


94 


16 


1.7 


94 zb 18 


32 


zb 


28 


4.5 


220576226 


32 zb 6 


283 zb 3 


1 


zb 


1 


30 zb 1 


30 


4 


1.1 


30 zb 12 


51 


zb 


14 


6.1 


220578040 


27 zb 9 


103 zb 2 




18.5 




654 


35 


16 


5.1 


247 zb 205 


50 


zb 


23 


7.0 


220584167 


59 zb 1 


178 zb 1 


9 


7 zb 


2 


234 zb 2 


24 


15 


7.9 


234 zb 35 


49 


zb 


21 


13.1 


220596913 


80 =b 2 


247 zb 1 





9 zb 


1 


141 zb 2 


153 


22 


4.3 


141 zb 10 


38 


zb 


28 


9.3 


910193711 


60 zb 8 


39 zb 4 


3 


5 zb 


7 


63 zb 8 


18 


17 


2.6 


63 zb 12 


80 


zb 


37 


4.1 


910279515 


52 zb 4 


108 zb 7 




3.8 




265 


70 


34 


5.3 


186 zb 14 


47 


zb 


40 


2.7 



The parameters of the kinematics modeling are the inclination (2), the position angle of the major axis (3), the turnover radius 
(4) and velocity (5). The inner slope is given in column (6). The mean residual of the velocity field is given in column (7) and the 
fit chi square in column (8). The maximum rotational velocity (9) and the mean velocity dispersion corrected from beam smearing 
(10) are computed after from the results of the fit. The extent of the velocity field (11) is derived using a S/N threshold of 3. 
The suffix "s" refers to secondary objects which have been detected in Ha and large enough to perform kinematics modeling. No 
error is given for r t and Vt when the plateau is not reached (see Section [3.2.3 ). 
1 Vmax is not coming from the model but from a detailed analysis (see Appendix [aI . 
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Table 4. Kinematics and close environment classification of MASSIV "first epoch" sample galaxies 



VVDS ID 


S/N 


^last 
Seeing 


APA [°] 


h e a r 1 2 




Shear 


Rotator 


Isolated 


Isolation flag 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


020106882 


7.6 


1.2 


23 ±4 


0.07 


3.2 


High 


Yes 


Yes 


B 


020116027 


8.2 


1.3 


23 ±6 


0.27 


0.6 


Low 


No 


No 


A 


020126402 


< 3 


- 


- 


- 


- 


- 


- 


- 


- 


020147106 


12.2 


1.4 


7±3 


0.18 


0.3 


Low 


Yes 


Yes 


B 


020149061 


6.2 


0.7 


34 ±90 


0.17 


1.5 


Low 


Yes 


Yes 


B 


020164388 


11.0 


1.2 


65 ±5 


0.11 


1.5 


Low 


No 


Yes 


B 


020167131 


3.7 


0.3 


87 ± 11 


0.22 


5.0 


Low 


No 


No 


A 


020182331 


6.1 


0.9 


34=1=6 


0.09 


2.0 


High 


No 


Yes 


B 


020193070 


4.9 


0.8 


31 ±3 


0.10 


3.6 


High 


No 


Yes 


B 


020208482 


4.1 


0.3 


20 ± 12 


0.18 


22.9 


Low 


Yes 


Yes 


B 


020214655 


7.5 


0.8 


51 ±5 


0.17 


0.8 


Low 


No 


Yes 


B 


020217890 


< 3 


- 


- 




- 


- 


- 


- 


- 


020239133 


5.3 


0.8 


14 ±6 


0.1 1 


2.0 


High 


Yes 


Yes 


B 


020240675 


5.5 


0.5 


7± 18 


0.40 


1.5 


Low 


No 


Yes 


B 


020255799 


4.5 


0.6 


75 ±77 


1.04 


0.2 


Low 


No 


Yes 


B 


020261328 


7.8 


1.0 


8±8 


0.07 


2.4 


High 


Yes 


Yes 


B 


020278667 


4.1 


0.2 


35 ± 14 


0.34 


1.5 


Low 


No 


Yes 


C 


020283083 


7.0 


0.9 


58 ±5 


0.18 


1.5 


Low 


No 


No 


B 


020283830 


4.8 


1.2 


14 ±2 


0.08 


11.3 


High 


Yes 


No 


B 


020294045 1 


7.0 


1.2 


2±4 


0.12 


3.9 


Low 


No 


No 


B 


020306817 


< 3 


- 


- 




- 


- 


- 


- 


- 


020363717 


11.4 


1.1 


49 ±8 


0.19 


0.5 


Low 


No 


Yes 


B 


020370467 


5.5 


0.9 


31 ±27 


0.45 


0.6 


Low 


No 


Yes 


B 


020386743 


8.8 


0.9 


67 ±8 


0.26 


0.8 


Low 


No 


No 


A 


020461235 


5.8 


1.1 


19 ±4 


0.10 


3.5 


High 


Yes 


No 


B 


020461893 


7.0 


1.3 


4±5 


0.10 


0.9 


Low 


Yes 


Yes 


B 


020465775 


6.8 


0.7 


31 ±5 


0.24 


0.8 


Low 


No 


No 


B 


140083410 


6.0 


1.0 


88 ± 12 


0.35 


0.4 


Low 


No 


Yes 


A 


140096645 


10.3 


1.0 


20 ±90 


0.07 


3.9 


High 


Yes 


No 


B 


140123568 


3.5 


0.1 


43 ±90 


0.44 


0.7 


Low 


No 


Yes 


B 


140137235 


3.2 


0.1 


8 ±29 


0.11 


3.5 


Low 


Yes 


Yes 


B 


140217425 


6.8 


1.9 


2± 1 


0.04 


7.1 


High 


Yes 


Yes 


B 


140258511 


8.1 


1.3 


52 ±42 


0.12 


5.1 


High 


Yes 


Yes 


A 


140262766 


6.5 


1.0 


32 ±28 


0.07 


3.1 


High 


Yes 


Yes 


B 


140545062 


8.5 


1.3 


14 ±41 


0.06 


3.1 


High 


Yes 


Yes 


B 


220014252 


11.0 


1.7 


5±2 


0.08 


1.4 


High 


Yes 


Yes 


B 


220015726 


10.4 


1.0 


8± 12 


0.04 


3.7 


High 


Yes 


Yes 


B 


220071601 


< 3 


- 


- 


- 


- 


— 


- 


- 


- 


220148046 


4.0 


0.4 


25 ± 18 


0.39 


0.9 


Low 


No 


Yes 


C 


220376206 


12.6 


2.4 


10 ± 1 


0.05 


2.8 


High 


Yes 


No 


B 


220386469 


4.5 


1.1 


17 ±12 


0.42 


0.9 


Low 


No 


Yes 


B 


220397579 


15.7 


2.0 


31 ±10 


0.68 


0.2 


Low 


No 


No 


A 


220397579s 


5.4 


1.2 


11±3 


0.06 


8.3 


High 


Yes 


No 


A 


220544103 


10.4 


1.2 


11 ±55 


0.07 


1.9 


High 


Yes 


No 


B 


220544394 


10.5 


1.1 


54 ±7 


0.08 


1.1 


Low 


No 


No 


A 


220576226 


11.6 


1.3 


45 ± 15 


0.15 


0.6 


Low 


No 


Yes 


B 


220578040 


5.9 


1.4 


10 ± 14 


0.08 


4.9 


High 


Yes 


Yes 


C 


220584167 


13.8 


2.1 


15 ±2 


0.04 


4.8 


High 


Yes 


Yes 


B 


220596913 


5.7 


6.2 


8± 1 


0.09 


3.7 


High 


Yes 


Yes 


B 


910193711 


6.0 


1.8 


30 ±4 


0.22 


0.8 


Low 


No 


Yes 


C 


910279515 


3.5 


1.5 


23 ±8 


0.11 


3.9 


High 


Yes 


Yes 


C 



Column (2) gives the average S/N of the S/N map. The size of the galaxy with respect to the seeing is given in column (3). The 
criteria used to determine if galaxies are in rotation or not (8) are the mismatch between morphological and kinematic position 
angles (4) and the residuals normalized by the velocity shear (5). The ratio of the rotation velocity over the local velocity dispersion 
is given in column (6). Galaxies are classified according the observed velocity shear in column (7): galaxies with V s hear < 100 
km s _1 or V s hear < 100 km s _1 are respectively classified as low and high shear). The close environment classification and its 
associated quality flag are given in (9) and (10). 

The suffix "s" refers to secondary objects which have been detected both in Ha and in the I-band image. 

1 This galaxy is classified as non-rotating even if it fulfill the criteria because the kinematics seem to indicate that this system is 
composed of two close companions (cf. Appendix |A| . 
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Table 5. Physical properties of MASSIV subclasses 



Class 


N 




star 






SFR 












a 










[in 10 M^l 






L ivi o y 1 


1 




[KpCj 






[km 








Med 




Dev 


Med 


IVlcdll 


Dev 


Med 




Dev 


Med 




Dev 


Rotating 


19 


3.9 


4.0 


3.7 


46 


60 


44 


3.1 


3.8 


2.4 


62 


58 


19 


Non- Rotating 


16 


1.2 


1.6 


1.2 


30 


39 


29 


2.7 


2.7 


1.2 


60 


61 


17 


Detected but unclassified 


10 


1.5 


2.3 


2.1 


24 


34 


30 


2.3 


2.3 


1.0 


43 


39 


24 


Isolated 


28 


1.4 


3.0 


3.4 


37 


47 


41 


2.7 


3.1 


2.2 


62 


56 


22 


Interacting 


13 


1.7 


2.1 


1.4 


30 


43 


32 


4.0 


3.9 


1.4 


53 


52 


21 


Rotating + Isolated 


14 


3.9 


4.1 


4.3 


58 


63 


49 


2.9 


3.8 


2.7 


62 


57 


19 


Rotating + Interacting 


4 


4.6 


3.6 


1.4 


64 


55 


33 


5.3 


4.2 


1.7 


73 


61 


25 


Non-Rotating + Isolated 


8 


1.4 


2.0 


1.7 


42 


37 


22 


1.9 


2.0 


1.1 


66 


64 


19 


Non-Rotating + Interacting 


7 


1.2 


1.3 


0.5 


26 


35 


33 


3.4 


3.5 


0.7 


53 


56 


15 



N is the number of galaxies in each subclass; med: median value; mean: average; dev: standard deviation. Only the main objects 
are considered. 
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Appendix B: Kinematic maps 

This appendix shows the kinematic maps for all detected galaxies. There are several sets of maps for galaxies with 
detected companions: one for the whole system (without kinematics modeling) and one for each modeled component. 
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Fig.B.l. Maps for VVDS020106882. From left to right: (top) the I-band CFHT image (arbitrary scale), the Ha flux 
map (arbitrary scale) and the S/N map, (middle) the observed velocity field, the rotating-disk-modeled velocity field, the 
residual velocity field, the uncertainty map on the velocity field, (bottom) the observed uncorrected velocity dispersion, 
the velocity dispersion map deduced from the velocity field model (beam-smearing effect and spectral PSF), the beam- 
smearing-corrected velocity dispersion map and the uncertainty on the velocity dispersion. The redshift is indicated in 
the top- left. In each map, north is up and east is left. The center used for kinematics modeling is indicated as a double 
black and white cross, the position angle is indicated by the black line. This line ends at the effective radius. The seeing 
FWHM is indicated on the Ha map as a circle. 
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Fig. B.2. Maps for VVDS020116027 system. From left to right: (top) the I-band CFHT image (arbitrary scale), the Ha 
flux map (arbitrary scale) and the S/N map, (middle) the observed velocity field, the uncertainty map on the velocity 
field, (bottom) the observed uncorrected velocity dispersion and the uncertainty on the velocity dispersion. The redshift 
is indicated in the top- left. In each map, north is up and east is left. The seeing FWHM is indicated on the Ha map as 
a circle. 
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Fig.B.3. Maps for VVDS020116027. Same caption as Figure B.l 
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Fig.B.4. Maps for VVDS020147106. Same caption as Figure B.l 
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Fig.B.5. Maps for WDS020149061. Same caption as Figure B.l 
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Fig. B.6. Maps for VVDS020164388. Same caption as Figure 



B.l 
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Fig. B.7. Maps for VVDS020167131 system. Same caption as Figure B.2 
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Fig.B.8. Maps for WDS020167131. Same caption as Figure 



B.l 
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Fig.B.9. Maps for WDS020182331. Same caption as Figure B.l 
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Fig. B. 11. Maps for VVDS020208482. Same caption as Figure 



B.l 
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Fig. B. 12. Maps for VVDS020214655. Same caption as Figure B.l 
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Fig. B. 13. Maps for VVDS020239133. Same caption as Figure B.l 
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Fig. B. 14. Maps for VVDS020240675 system. Same caption as Figure 



B.2 
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Fig. B. 15. Maps for VVDS020240675. Same caption as Figure B.l 
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Fig. B. 16. Maps for VVDS020255799. Same caption as Figure B.l 



j 



- SNR 














VVDS020261328 
z = 1.5290 




Aa ("] 





Aa (-; 



V r residuals 



J_ 




Fig. B. 17. Maps for VVDS020261328. Same caption as Figure B.l 



■J 



20 _! - 




B. Epinat et al.: MASSIV: Mass Assembly Survey with SINFONI in VVDS 



VVDS020278667 
z = 1.0516 




Fig. B. 18. Maps for VVDS020278667. Same caption as Figure B.l 
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Fig. B. 19. Maps for VVDS020283083. Same caption as Figure B.l 
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Fig. B.20. Maps for VVDS020283830 system. Same caption as Figure B.2 
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Fig. B. 21. Maps for VVDS020283830. Same caption as Figure B.l 
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Fig. B. 22. Maps for VVDS020294045. Same caption as Figure B.l 
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Fig. B. 23. Maps for VVDS020363717. Same caption as Figure B.l 




Fig. B.25. Maps for VVDS020386743 system. Same caption as Figure 



B.2 
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Fig. B. 26. Maps for VVDS020386743. Same caption as Figure B.l 
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Fig. B. 27. Maps for VVDS020461235. Same caption as Figure B.l 
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Fig. B. 28. Maps for VVDS020461893. Same caption as Figure B.l 
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Fig. B. 30. Maps for VVDS140083410. Same caption as Figure B.l 
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Fig. B. 31. Maps for VVDS140096645 system. Same caption as Figure B.2 






Fig. B.33. Maps for VVDS140123568. Same caption as Figure 



B.l 
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Fig. B.34. Maps for VVDS140137235. Same caption as Figure 
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Fig. B. 35. Maps for VVDS140217425. Same caption as Figure B.l 
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Fig. B. 36. Maps for VVDS140258511. Same caption as Figure B.l 
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Fig. B. 37. Maps for VVDS140262766. Same caption as Figure B.l 
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Fig. B. 38. Maps for VVDS140545062. Same caption as Figure B.l 
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Fig. B. 39. Maps for VVDS220014252. Same caption as Figure B.l 



44 



B. Epinat et al: MASSIV: Mass Assembly Survey with SINFONI in VVDS 




Fig.B.41. Maps for VVDS220148046. Same caption as Figure 



B.l 
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Fig. B.43. Maps for VVDS220376206. Same caption as Figure 



B.l 




Fig. B.45. Maps for VVDS220397579 system. Same caption as Figure 



B.2 




Fig. B.47. Maps for VVDS220397579 companion. Same caption as Figure 



B.l 
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Fig.B.48. Maps for VVDS220544103. Same caption as Figure B.l 
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Fig. B.49. Maps for VVDS220544394 system. Same caption as Figure B.2 
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Fig. B. 50. Maps for VVDS220544394. Same caption as Figure B.l 
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Fig. B. 51. Maps for VVDS220576226. Same caption as Figure 



B.l 
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Fig. B.52. Maps for VVDS220578040. Same caption as Figure 
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Fig. B. 53. Maps for VVDS220584167. Same caption as Figure B.l 
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Fig. B. 54. Maps for VVDS220596913. Same caption as Figure B.l 
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Fig. B. 55. Maps for VVDS910193711. Same caption as Figure B.l 
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Fig. B.56. Maps for VVDS910279515. Same caption as Figure 



B.l 
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